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Abstract 



Recent inapproximability results of Sly (2010), together with an approximation algorithm 
presented by Weitz (2006) establish a beautiful picture for the computational complexity of 
approximating the partition function of the hard-core model. Let Ac(Ta) denote the critical 
activity for the hard-model on the infinite A-regular tree. Weitz presented an FPTAS for the 
partition function when A < Ac(Ta) for graphs with constant maximum degree A. In contrast, 
Sly showed that for all A > 3, there exists ea > such that (unless RP = NP) there is no 
FPRAS for approximating the partition function on graphs of maximum degree A for activities 
A satisfying Ac(Ta) < A < Ac(Ta) -I- ea- 

We prove that a similar phenomenon holds for the antiferromagnetic Ising model. Sinclair et 
al. recently extended Weitz's approach to the antiferromagnetic Ising model, yielding an FPTAS 
for the partition function for all graphs of constant maximum degree A when the parameters of 
\^ , the model lie in the uniqueness region of the infinite A-regular tree. We prove the complementary 

P\J ' result that for the antifcrrogmanctic Ising model without external field that, unless RP = NP, 

Cn , for all A > 3, there is no FPRAS for approximating the partition function on graphs of maximum 

degree A when the inverse temperature lies in the non- uniqueness region of the infinite tree Ta . 
CO . Our proof works by relating certain second moment calculations for random A-regular bipartite 

^^ ' graphs to the tree recursions used to establish the critical points on the infinite tree. 
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1 Introduction 

A remarkable computational transition has recently been established for the complexity of approx- 
imating the partition function of the hard-core model. Amazingly, this computational transition 
coincides exactly with the statistical physics phase transition on infinite A-regular trees. In this 
work we prove that a similar phenomenon holds for a much more general class of 2-spin models. 

Our general setup is 2-spin systems on an input graph G = {V, E) with maximum degree A. 
We follow the setup of several related previous works [10, 23, 16]. Configurations of the system are 
assignments a : V ^ {— 1,-|-1}. For a given configuration a G {— 1,-|-1}^, denote by n-{a) the 
number of vertices assigned —1, m~{a) the number of edges with both endpoints assigned —1 and 
m~^{a) the number of edges with both endpoints assigned -|-1. 

There are three non- negative parameters Bi, B2 and A where Bi is the edge activity for (— , — ) 
edges, B2 is the edge activity for (-|-,-|-) edges, and A is the vertex activity (or external field). 
Configuration a G { — 1,-|-1} has weight 



u;g(c7) = A"-MBi"^"(")52"^^("\ 



The Gibbs distribution // is over the set {—1, +1}^ where ij{a) = WG{cr)/Z, where Z = Zg{Bi,B2, A) 
is a normalizing factor, known as the partition function, defined as: 

z:= Yl yjGi'j)= Yl A"-Mi?r"("^i?2"^("^. 

(Te{-i,+i}^' o-e{-i,+i}^' 

Let us point out a few examples of 2-spin models that are of particular interest. Setting 
Bi = and B2 = 1, the only configurations with positive weight in the Gibbs distribution induce 
independent sets of G. This case is known as the hard-core model with fugacity A. Setting Bi = 
B2 = B, one recovers the Ising model. The case A = 1 is the Ising model without external field. 
When i? < 1 it is the antiferromagnetic Ising model, and when i? > 1 it is the ferromagnetic Ising 
model. In general, when B1B2 < 1 the model is called antiferromagnetic, and when B1B2 > 1 it is 
ferromagnetic. 

The Ising and hard-core models (and more generally spin models) have been studied thoroughly 
in various contexts, perhaps more persistently in Statistical Physics. Their interpretation as ide- 
alized models of microscopic interaction within a body have for a long time nurtured research on 
understanding their macroscopic behavior on appropriate underlying graph structures such as the 
infinite grid or the infinite triangular lattice. The intriguing questions arising have made an impact 
on other fields as well, most notably Probability and Computer Science. For illuminating accounts 
of such directions, we refer the reader to [18] and [26]. For the sole purpose of illustrating a toy 
relevant example, a moment's thought reveals that even determining a ground state of the Gibbs 
distribution, i.e., a configuration with largest weight, reduces to NP-complete problems in the cases 
of the Ising antiferromagnetic and hard-core models. 

In this paper, we focus on the computational complexity of approximating the partition function 
in 2-spin models. A long series of works, which we only touch upon later in this introduction, have 
identified the connection of the problem with decay of correlation properties in Gibbs measures of 
the infinite A-regular tree. In statistical physics terms, this is known as the phase transition on 
the Bethe lattice. To briefly sketch this important concept, let Ta denote the infinite A-regular 
tree. When the parameters of the 2-spin model lie in the so-called uniqueness region of Ta, for 
finite complete trees of height h, the influence on the root from fixing a configuration on the leaves 
dies off in the limit as the height goes to infinity. In sharp contrast, in the non-uniqueness region 



there exist configurations of the leaves for which their influence on the root persists in the limit. 
For convenience, we often refer to the uniqueness region omitting the reference to Ta- 

The uniqueness/non-uniqueness regions can be determined rather easily on Ta (for other infinite 
graphs this is far from easy, see [3] and [22] for recent advances in the two dimensional grid). For 
the hard-core model, Kelly [15] showed non-uniqueness holds iff A > Ac(Ta) := (a-2)^ ' ^'~'^ ^^^ 
antiferromagnetic Ising model without external field, non-uniqueness holds iff i? < Bc{T/\) := 
-^- (see, e.g., [23]). For general antiferromagnetic 2-spin models with soft constraints (i.e., 
B1B2 > 0), non-uniqueness holds iff y/BiB2 < (A — 2)/A and A G (Ai,A2) for some critical 
values Ai(A,Si,52),A2(A,5i,S2) (see, e.g., [16]). 

Before stating our results, it will be instructive to review some known results from the literature. 
When A = 1, the partition function is the number of independent sets in G, and this quantity is 
T^P-complete to compute exactly [27], even when A = 3 [11]. For the Ising model, once again, 
exact computation of the partition function for graphs of maximum degree A is #P-complete, 
even for the ferromagnetic model [14]. Hence, the focus is on the computational complexity of 
approximating the partition function. 

This is where the phase transition on Ta comes into play. Specifically, for the hard-core model 
it was not until recently that a deep connection between uniqueness/non-uniqueness and the com- 
plexity of approximating the partition function was established. Namely, for graphs of constant 
maximum degree A, Weitz [28] presented a beautiful FPTAS for approximating the partition func- 
tion when A < Ac(Ta). On the other hand, Sly [24] presented a clever reduction which proves that 
for every A > 3, there exists sa > 0, such that it is NP-hard (unless RP=NP) to approximate the 
partition function for graphs of maximum degree A for any A where Ac(Ta) < A < Ac(Ta) -|- sa- 

It was believed that Sly's inapproximability result should hold for all A > Ac(Ta). However, 
Sly's work utilized results of Mossel, Weitz and Wormald [21] which showed that for a random 
A-regular bipartite graph, an independent set chosen from the Gibbs distribution is "unbalanced" 
with high probability for A under the same condition as in Sly's result. 

Mossel et al.'s work [21] involves a technically complicated second moment argument that 
Sly [24] calls a "technical tour-de- force" . In [7] we extended Mossel et al.'s work to all A > Ac(Ta) 
for A = 3 and A > 6. However, that work built upon Mossel et al.'s work and had an even more 
complicated analysis. In this work, as a byproduct of the approach we devise for attacking the Ising 
model, we can prove inapproximability for the hard-core model for all A > Ac(Ta) for A = 3,4,5, 
and our proof is simpler than in [7]. This settles the picture for the hard-core model, and also 
proves Conjecture 1.2 of Mossel et al. [21]. 

Theorem 1. Unless NP = RP, for the hard-core model, for all A > 3, for all A in the non- 
uniqueness region of the infinite tree Ta, there does not exist an FPRAS for the partition function 
at activity A for graphs of maximum degree at most A. 

Our main focus in this paper is to address whether the above phenomenon for the hard-core 
model occurs for other models. In particular, our goal is to address for 2-spin antiferromagnetic 
models whether the computational complexity of approximating the partition function for graphs 
with maximum degree A undergoes a transition that coincides exactly with the uniqueness/non- 
uniqueness phase transition on the infinite tree Ta. 

Let us first review existing results in the literature. For the ferromagnetic Ising model, Jerrum 
and Sinclair [14] presented an FPRAS for all graphs, for all -B > 1 where B = Bi = B2. For 
the antiferromagnetic Ising model, for general graphs there does not exist an FPRAS for B < \ 
[14], unless NP = RP. For constant A, Sinclair et al. [23] extend Weitz's approach to the 
antiferromagnetic Ising model (see also, [30, 22, 16] for other results on the Ising model), yielding 
in the uniqueness region an FPTAS for the partition function of graphs with maximum degree A. 



We prove for the antiferromagnetic Ising model without external field that for all B in the 
non- uniqueness region of the infinite tree there does not exist an FPRAS for the antiferromagnetic 
Ising model. 

Theorem 2. Unless NP = RP, for the antiferromagnetic Ising model without external field, for 
all A > 3, for all B in the non-uniqueness region of the infinite tree T^, there does not exist 
an FPRAS for the partition function at inverse temperature B for graphs of maximum degree at 
most A. 

The hard-core model and Ising model are the two most well-studied examples of 2-spin systems. 
For ferromagnetic 2-spin systems with no external field (i.e., B1B2 > 1 and A = 1), Goldberg et 
al. [10] presented an FPRAS for the partition function for any graph. For any antiferromagnetic 
2-spin model, for constant A, an FPTAS for the partition function was obtained by Sinclair et al. 
[23] for A-regular graphs in the uniqueness region of the infinite tree Ta, and by Li et al. [16] for 
graphs of maximum degree A in the intersection of the uniqueness regions for infinite trees T^ for 
d < A. For antiferromagnetic 2-spin models we obtain a complementary inapproximability result 
in the non-uniqueness region. Our results do not reach the uniqueness/non-uniqueness threshold 
in general. 

Theorem 3. Unless NP = RP, for all d > 2, for all Bi,B2,X that lie in the non-uniqueness 
region of the infinite tree T^+i and \fB\B2 > rf , there does not exist an FPRAS for the partition 
function at parameters Bi,B2, X for graphs of maximum degree at most A = d + 1 . 

1.1 Independent Results of Sly and Sun 

In an independent and simultaneous work. Sly and Sun [25] obtained closely related results. They 
prove inapproximability of the partition function in the case of the hard-core model and the Ising 
model. Their result for the Ising model also covers the case of an external field (which then extends 
to general 2-spin models [23, 25]). 

The main technical result in our proof is a second moment argument to analyze the partition 
function on random A-regular bipartite graphs, as outlined in the following Section 2. As a con- 
sequence we can estimate the partition function within any arbitrarily small polynomial factor. In 
contrast. Sly and Sun's approach builds upon the interpolation scheme of Dembo et al. [5] to ana- 
lyze the logarithm of the partition function for random A-regular bipartite graphs from which they 
only obtain estimates within an arbitrarily small exponential factor. Since their proof yields weaker 
results regarding the partition function for random bipartite graphs, their NP-hardness reduction 
requires some modification of the approach appearing in [24] . 

Finally, we make some minor relevant remarks. Sly and Sun observe that their proof yields 
inapproximability of the partition function within an exponential factor. This also holds for the 
proof in [24] and hence in our case as well. 

2 Proof Outline 

The main element of the proof of Theorem 3 (and similarly, the proofs of Theorems 1 and 2) is 
to analyze random A-regular bipartite graphs to show a certain bimodality as in Mossel et al.'s 
work [21] for the hard-core model. From there the same reduction of Sly [24] applies with some 
non-trivial modifications in the analysis. Hence, in some sense our main technical result is to show 
that in random A-regular bipartite graphs, under the conditions of Theorem 3, a configuration 
selected from the Gibbs distribution is "unbalanced" with high probability. 



To formally state our results, let us introduce some notation. Let G{n, A) denote the probability 
distribution over bipartite graphs with n+n vertices formed by taking the union of A random perfect 
matchings. We will denote the two sides of the bipartition of the graphs as V\^V2- Strictly speaking, 
this distribution is over bipartite multi- graphs. However, since our results hold asymptotically 
almost surely (a.a.s.) over ^(n. A), as noted in [21], by contiguity arguments (see the note after 
the proof of [19, Theorem 4]), our results also hold a.a.s. for the uniform distribution over bipartite 
A-regular graphs. For a complete account of contiguity, we refer the reader to [13, Chapter 9]. 

For a graph G ~ ^(n. A) and a given configuration a G { — 1, +1}^, denote by N-{a) the set of 
vertices assigned —1 and by N^{a) the set of vertices assigned +1. For a,/3 > 0, let 



■^a,/3 



{cjE{-1,+1}^| \N_{a) nVil =an, | iV_ (a) n T/2 1 = /3n}, 



that is, a (resp. /3) is the fraction of the vertices in Vi (resp. V2) whose spin is —1. Denote by 
•£Bai |-|g ^j^g gg|. q£ ^Yl balanced configurations, and for 6 > 0, denote by T,^ the set of (5n-unbalanced 
configurations. Formally, 

a \a~l3\>5 

Theorem 4. Under the hypotheses of Theorem 3, there exist a > \ and 5 > 0, such that asymp- 
totically almost surely, for a graph G chosen from ^(n. A), the Gihhs distribution /i satisfies: 

Therefore, the Glauber dynamics is torpidly mixing. 

Theorem 4, while interesting in its own right, provides us also with the necessary ingredients 
to carry out Sly's reduction. To prove Theorem 4, we use a second moment approach as in [21]. 
Specifically, for G ~ Q{n,/S.) we analyze the random variable ^g , where Zq = 'Ylia&T.^'P ^g(c")) 
for some well-suited values of a, /3. Intuitively, the a, /3 should be selected so as to maximize the 
first moment of Zq . To have some hope in proving Theorem 4, it better be the case that the a, j3 
are not equal. This is indeed true in the non-uniqueness region and moreover there is an explicit 
connection with the extremal measures in the infinite tree. To formally depict this connection, let 
us make a short detour on the extremal measures in the infinite A-regular tree Ta. For a more 
complete account of extremality, we refer the reader to [9] . 

For every i?i,i?2,A > there exists a unique translation invariant Gibbs measure /i on Ta, 
known as the free measure with occupation density p* = iJ,{ap) for p the root of the tree. In the 
non- uniqueness region, there also exist two semi-translation invariant measures /i+ and /i_ whose 
occupation densities we denote by p~^ = iij^{ap), p~ = n^{ap). These measures are obtained by 
conditioning on level 2i (resp. 2i + 1) of the tree to be completely occupied and taking the weak 
limit as £ — )• 00. In our context, the phase transition on the infinite (A — l)-ary tree Ta rooted 
at p will also be of interest. The uniqueness/non-uniqueness regions for this slightly modified 
tree coincide with those of Ta, albeit with different occupation probabilities for the root in the 
analogous free and extremal measures, which we denote by q*,q^- Using monotonicity arguments, 
it can easily be proved that p~ < p* < p"*" and q~ < q* < q'^. 

Lemma 1. For the distribution Q{n,A), lim —\ogKg[Zp] is maximized for: 

n— >oo n 

1. (a,/3) = {p*,p*), whenever Bi,B2,X are in the uniqueness region of the infinite tree. 

2. (a,/3) = (p^,p^), whenever Bi, B2, X are in the non-uniqueness region of the infinite tree. 



Note that Lemma 1 already yields a weak form of a bimodality in the Gibbs distribution 
on random A-regular graphs: in the non-uniqueness region, the configurations with the largest 
contribution are unbalanced in expectation. Since Zq is typically exponential in n, this already 
provides a high probability upper bound on the number of balanced configurations (using Markov's 
inequality) . 

However, to get tight results, we need a strong form of this bimodality. To do this, we look more 
carefully at the random variable Zq for (a, /3) = (p^,p^). Note that if this variable is concentrated 
around its expectation, Lemma 1 establishes that it dominates exponentially variables Z^ with 
||(a',/3') — (p ,p^)|| > e where e can be made arbitrarily small for sufficiently large n. Hence, one 
obtains the required bottleneck that is needed for the proof of Theorem 4. The torpid mixing of 
Glauber dynamics follows then by a conductance argument given in [6]. 

In fact, the same idea yields something stronger: a typical configuration of the system will 
be from the set YP '^ (roughly). Informally, this allows us to say that with high probability G 
behaves as a boolean gadget depending on the phase of a configuration, where the phase is the 
bipartition which has the largest number of vertices assigned -1. Pasting a modified version of G 
(while maintaining a variant of this boolean behaviour) on the vertices of a graph if, it turns out 
that the partition function of H is dominated by configurations where the phases of the gadgets 
correspond to a MaxCut of H. Of course, the quantification of this nice idea requires serious 
effort, but this has already be done in [24]. To hint why such a scheme succeeds, note for example 
that in the antiferromagnetic Ising model, a ground state of the Gibbs distribution corresponds 
to a MaxCut. Thus, the reduction ensures that the phase configurations induced by MaxCut 
partitions of the original graph K dominate exponentially the partition function. 

Let us now return to the technical core of the method, which requires analyzing the second 
moment of Z^ for (a,/?) = (p ,p^)- By symmetry, we may clearly focus on (a,/3) = (p''',p~)- Pick 
two configurations cri,(T2 from the set TP 'P (not necessarily distinct). The two most important 
variables here are the overlap of the configurations o"i , o"2 . Formally, let 

Ks = {(^1' ^2)|cTi,a2 G s"'^ |iV_((Ti) n iv_(a2) nVi\=-f, |iv„((7i) n iv„(cj2) n V2\ = 6} , 

and yj' = J2( \py;a,i3 wg{o'i)wg{(J2)- We will study yj' for {a, (3) = {p~^,p~), so this is why 

we dropped the dependence on a, f3. Clearly, Eg[Y'J' ] is the contribution of pairs of configurations 
in Ti'^'P with overlap "fn,6n to the second moment Eg[(Z^ )^]. 

In our case, the second moment approach has any chance of succeeding when the largest con- 
tribution to the second moment comes from uncorrelated pairs of configurations. This is captured 
by the following condition. 

1 5 

Condition 1. For (a, /3) = (p"'~,p~), lim —logKglYp,' ] is maximized at 'f = 0^,6 = /3^. 

n— >oo n 

Lemma 2. Under the hypotheses of Theorem, 3, Condition 1 holds. 

Lemma 3. For the antiferromagnetic Ising model with A = 1, A = 3 and B in the non-uniqueness 
region, Condition 1 holds. 

Lemma 4. For the hard-core model, A = 3, 4, 5 and A in the non-uniqueness region. Condition 1 
holds. 

The innocuous statements of Lemmas 2, 3, and 4 are rather misleading. Namely proving that 
uncorrelated pairs of random variables "dominate" the second moment is a standard complication 



of the second moment approach in a variety of settings, see for example [1, 2, 21]. To highhght the 
technical difficulties, the second moment typically has a number of extra other variables that need 
to be considered and those are usually treated by sophisticated analysis arguments. We escape this 
paradigm in the following sense: our arguments try to explain the success of the second moment 
by relating it to the intrinsic nature of the problem at hand (in our context, uniqueness on the 
infinite tree). To stress the gains of such a scheme, it underlies perhaps a more generic method to 
translate uniqueness proofs to similar second moment arguments. 

In the hard-core model, due to the hard constraints, only one extra variable (besides 7, 5) was 
needed and it could be expressed explicitly in terms of 7, 6, so that, in fact, one had a 2-variable 
optimization. Even with this seemingly simple setup and despite the variety and complexity of 
approaches appearing in [21], [24] and [7], certain regions for A = 4, 5 remained open. In the 
general 2-spin model the second moment is more complex than in the hard-core model in the sense 
that the soft constraints cause the number of extra variables to blow-up into nine. To make matters 
even worse, one cannot get rid of those extra variables in some easy algebraic way, since such a 
scheme reduces to a sixth order degree equation in terms of 7, 5. These obstacles are surmounted 
in the proof of Lemma 2 by a very clean argument: the analysis of the second moment is reduced 
to a certain stable form of tree recursions on Ta, allowing for a short inequality argument. 

After establishing Condition 1, one obtains that for (a, /3) = (p , p^) it holds that Eg[(Z^ )^] ~ 

^g[^g ]• ^^ ^^^^ point, we are in position to study the ratio Eg[(Z^ )^]/(Eg[Z^ ])^. This is 
a technical computation and requires an asymptotic tight approximation of sums by appropriate 
Gaussian integrals. While the asymptotic value of the ratio is a constant, as in [21], the constant 
is greater than 1. 



Lemma 5. When Condition 1 is true, for (a,/3) = {p ^p^), it holds that 

Eg 



lim 



[Z^d'f 



where uj is given by (4). 



{i-uy^^~'^^\i-{A-i)Wy'/\ 



Remark 1. The formula given in the lemma is valid for the hard-core model as well (Bi = 0, i?2 = 
1), agreeing with the statement of [21, Theorem 3.3]. 

This complication does not allow to get asymptotically almost surely results using the second 
moment method. This is by now a standard scenario in the theory of random graphs and it can be 
soothed by the small graph conditioning method [13, 29]. The application of the method utilizes 
the moments' analysis. In our setting, the method explains the variance of the variables Zq 
by conditioning on the occurrence of a cycle in a random A-regular graph. The consequence of 
applying the small graph conditioning method is the following result. 

Lemma 6. When Condition 1 is true, for (a,/3) = {p^^p"^), it holds asymptotically almost surely 
over the graph G ~ Q{n, A) that 

Using Lemmata 1 and 6, one can carry out the proof of Theorem 4 as in [21, Proof of Theorem 
2.2] (see also [7, Proofs of Theorems 2 and 3]). This remaining argument was sketched earlier in 
this section. 

The outline of the rest of the paper is as follows. In Section 3 we derive certain properties 
for the infinite tree which are used throughout the remainder of the paper. We then look at the 

6 



first and second moments in Section 4. In Section 5 we look at the optimization problems for 
the first and second moments and prove Lemmata 1, 2. In Section 6 we outline the proof of our 
inapproximability results based on Sly's reduction for the hard-core model. 

In Section 7 we analyze the ratio of the second moment to the first moment squared for random 
A-regular graphs and for the modified graph used in Sly's reduction. We then apply the small 
graph conditioning method in Section 8 to get asymptotically almost surely results. Finally, in 
Section 9 we give the proofs of some remaining lemmas, including Lemmata 3 and 4. 

3 Extremal Measures on the Infinite Tree for 2- Spin Models 

In Section 2, we defined the measures on T^ and Ta which will be of interest to us. Here we state 
some properties for the densities q^ , q* and p^,p*, which we are going to utilize throughout the 
text. 

Standard tree recursions for Gibbs measures (e.g., see [23]) establish that 



1^ _x/^^-T^ + M^' 9* _x^^l-T^ + l^''"' 



?^ \ j^+B2 / ' 1-9* V T^ + ^2 



(1) 



,± 



It will be very handy to have separate expressions for j^-^, j^r^, i-e., 

1 — q+ 1 — ? I — q* 

The following lemma summarizes the properties of the occupation densities g"*", q^ ,q* which we are 
going to utilize. Its rather folklore proof is a straightforward analysis of appropriate tree recursions 
and is omitted. 

Lemma 7. For A > 3 and Bi,B2,X > 0, consider the system of equations 

with x,y > 0. 

1. In the uniqueness region, the only solution to (3) is {Q*,Q*). 

2. In the non-uniqueness region, (3) has exactly three solutions which are exactly {Q^,Q^) and 
{Q*,Q*). It holds that Q~ < Q* < Q+ . 

We also define the following quantities w, uj* which will emerge throughout the paper. 

^ (1 - BiB2fQ+Q- , ^ il-BiB2)HQ*? ... 

^' iB2 + Q-){B2 + Q+)il + BiQ-)il + B2Q+)' '^' {B2 + Q*f{l + BiQ*Y' ^' 

Note that w* is obtained by w after identifying Q with Q* . The following lemma for uj,uj* will 
be crucial to prove some optimality conditions required to establish Lemma 1 and Condition 1. Its 
proof is given in Section 9.1. 

Lemma 8. For {Bi,B2,X) in the non-uniqueness region, it holds that 

(A-l)2a;<l<(A-l)V. 

7 



Let us explain the mysterious looking inequality in Lemma 8. While not obvious, it is related 
with the two-step recurrence on the infinite tree Ta- The two step-recurrence is obtained by the 
composition of two successive one step recurrences. In the non-uniqueness region, the fixed points 
of the two step recurrence are exactly q^ , q* (this is basically Lemma 7) . The inequality in Lemma 8 
establishes, after some calculations, that the fixed point q~^ (resp. q~) is an attractor of the two step 
recurrence in /i"'" (resp. fi,~), while the fixed point q* is repulsive. More specifically, the derivative 
of the two step function at q is equal to (A — l)^a;, and similarly at q* is equal to (A — l)^a;*. 

This establishes at first place that the infinite volume measures /x"*" , //" are well defined (through 
the DLR condition, see [9]) and moreover it can be utilized to prove a certain form of non- 
reconstruction, see for instance [20], in the measures fi . From a measure-theoretical viewpoint 
these facts come down to the following: in the non-uniqueness region, the infinite volume measures 
on the infinite tree form a simplex with vertices fi . We should also note that the analysis in 
[23] (see also [16]) establishes the following criterion for uniqueness, which we reformulate for our 
purposes: uniqueness holds on Ta iff (A — l)^a;* < 1. 

Finally, let us translate the above to obtain a handle on p ,p*. Using a slight modification of 
the tree recursions in Ta, one obtains that 



p± _ g± Pl-T^ + 1\ p* _ q* B,.j^ + 1 



l_p± 1 - g± \ _il_ + 5 I ' 1 - P* l-q* \ j3^ + B2 

Using (2) and (5), an easy calculation gives 



(5) 



Equation (6) will prove to be very useful to rationalize and simplify seemingly intricate calculations. 
For the case of the Ising model with zero external field {Bi = B2 = B,X = 1), it holds that 
p'^ +p~ = q~^ + q~ = 1, p* = q* = 1/2. These properties may be proved easily using the symmetries 
of the model. Namely, observe that if {x,y) satisfy (3) then so do (1/y, I/2;). By items 1 and 2 of 
Lemma 7, this implies that q^ + q~ = 1 and q* = 1/2. By (5), it also follows that p~^ + p~ = 1, 
p* = 1/2. 

4 Moment Analysis 

In this section, we give expressions for Kg [-^q ] , Eg [Yq ] and then give the asymptotics of 



^logE,[Z^'^],ilogE,[yJ'^l 



an\ /(I— o)n^ 



4.1 First Moment 

The first moment of Zq is given by 

Eg [Zf] = A"("+'^) ('')(!!]( Y. ^^"^^|f-^^^"^ i?f»i?(^-°-^+^)" ) . (7) 

\anj \pnj \^5<„ ,5<^ (^J J 

Ignoring the contribution of the external field A, the terms outside the sum count |Sg | . The sum 
gives the expected contribution of o" G ^g for one matching, where Sn is the number of edges 
between N-{a) D Vi and N-{a) n V2 in a single matching. The weight of a for a single matching 



depends only on 6n and is equal to B1B2 " . The fraction in the sum is the probability 

that there are exactly 6n edges from N^{a) n V2 to A^4.(cr) n V2. Since the matchings are chosen 
independently, the exponentiation to the power A gives the expected distribution of o" S H^ for 
A matchings. Finally, the A term in the beginning accounts for the contribution of the external 
field for a G ^g ■ 

A different way to arrive at a slightly different formulation for the first moment follows. Here, we 
introduce some redundant variables at the cost of introducing some complexity of the presentation, 
but underlying a general approach to treat such expressions. 

Observe that the sets N^ {a)r]Vi, N~^ (a) CiVi form a partition of Vi and similarly for V2- Thus 
for a perfect matching let Xij, (1 < i,j < 2) denote the number of edges having their endpoints in 
partition i of Vi and in partition j of ¥2- Note that the Xij must satisfy 

xu+xi2 = a, xii+X2i=/3, 

aJ2i + X22 = 1 - a, X12 + X22 = 1 - /3, (8) 

Xij > 0. 



The first moment may be written as 



an \ f{l—a)n\ ( fin \ /{I— /3)n\ ^ ^ 



/ r, \ / r, \ f ( <^^) ('U-a]"^ ( P" ^ Al-/^)"^ 

Ig Iz'^^t^] = \ {"+/3)'" ( 111! \^ \xiinj V X2in J \xiinj V X2in J rjxnn T^X22n 1 /g\ 



^2 
\xiin,xi2n,X2in,X22n/ 



J 



In the above sum, we denote by X the Xij which satisfy (8). Note that the expressions (7) and 
(9) are completely equivalent after identifying 5 with xu and using (8) to express the variables 
2^12 5 2;2iia;22 in terms of xn. 

4.2 Second Moment 

We continue with the ideas presented in Section 4.1, to get an expression for Eg \Y^' 1 . Now there 
are four induced parts by cri, (T2 on each Vi and V2. Let 'jn, 6n stand for the number of vertices in 
Vi and V2 whose spin is —1 in both ai, a"2 and t/ijU the number of edges incident with endpoints in 
partition i of Vi and partition j of V2 in a single matching. For ease of exposition, define 

Li = 7, L2 = -L3 = a - 7, L4 = l-2a + 7, 
R^ = 5, R2 = Rs = (3 - 6, R4 = l-2l3 + 6. 

In the same spirit with expression (9) for the first moment, we have that 

Mrr] - A*-."(:„) (;) (;•;) (;-_;';,;) (t) (<-_«;) m 

Y 1=1 Vyii"-' Vi'i'n; Visn, ViAnJ ^^^ \xijn, X2jn, x^jU, x^juj \xiin, x^n, ..., x^in 

. A 

^ o(2j/ll+3/12+S/13+J/21+3/22+J/31+y33)»^ 73(S/22+3/24+J/33+?/34+S/42+J/43+2j/44)n. 

In the above sum, we denote by Y the yij which satisfy 



T,jyij=Li, fori = l,...,4 

Eiyij=Rj, forj = l,...,4 (11) 

Vij > 0. 
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The interpretation of formula (10) should be clear: ignoring the external field, the first line accounts 
for |S"'^|, a term in the sum accounts for the weight of a matching specified by the cardinalities Y 
(last line) as well as the probability of such a matching occurring (middle line) . 

4.3 Asymptotics of the Logarithms of the Moments 

We next turn to finding the asymptotics of -logEg[Z^ ], -log Eg [Y^' ]. The asymptotic order 
of the sums in (9) and (10) are dominated by their maximum terms. For the sole purpose of having 
a unifying treatment for soft and hard constrained 2-spin models, here and in the rest of the paper, 
we adopt the usual conventions that InO = — oo and OlnO = 0. Standard application of Stirling's 
approximation yields for the first moment: 

lim -logEcTZ^'^l =max$i(a,/3,X), (12) 



where <I>i(a,/3,X) 

/i(a,/3) 

91 (X) 



(a + /?) In A + (A - l)/i(a, /?) + A<7i(X) 

alna + (1 - a) ln(l - a) + /31n/3 + (1 - /3) ln(l - ^) 

(xii In^i + X22 lnB2) - Y^ij^ij Inxjj. 



And for the second moment: 



lim -logEcfy^^' 1 =max$2(7,<J,Y), (13) 

where ^2(7, S,Y) := 2{a + 13) In X + {A - 1) f 2(7 , 6) + Ag2{Y) 

/2(7, 6) := 2{a - 7) ln(a - 7) + 7 ln7 + (1 - 2a + 7) ln(l - 2a + 7)+ 
2(/3 - 6) ln(/5 - 6) + 5ln6+ (1-2^ + 6) ln(l -2/3 + 6) 

52 (Y) := (2yii + yu + yn + 2/21 + ^22 + ysi + yss) In Si 

+ (y22 + y24 + yss + y-u + y42 + y^s + '^yu) in B2 - J2i,jyij in yij. 

Recall that we are going to study the second moment for specific values of q, /3, namely (a, (3) = 
(p^,p^), so this is why we dropped the dependence of the functions $2, /2 on a, j3. The functions 
31,(72 are defined on the regions (8), (11) respectively. For fixed a,/3,7,(5, they are strictly concave 
over the convex regions on which they are defined and hence they have a unique maximum, so this 
also the case for the functions $1, $2- The limits (12) and (13) can thus be justified using standard 
Laplace arguments (see for example [4]). 

5 Finding the Maxima 

In this section, we use the information obtained in Section 4.3 to establish Lemmata 1, 2. The 
proofs for Lemmata 3 and 4 are given in Sections 9.5 and 9.6 respectively, see also Section 5.4 for 
an overview. These proofs are based on the approach given in this section. 

5.1 Optimizing Entropy Distributions 

By the limits (12) and (13), to establish Lemmata 1, 2 (and 3, 4 as well), it suffices to study the 
maxima of the functions 

4>i{a,j3) := max$i(a,/3,X) and(/)2(7,5) := max$2(75 (5, Y). 
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A small abstraction allows us to treat the two functions in a unifying way. Namely, observe 
that (j)ii4'2 have the form (j){u,v) = maxz $(n, -y, Z), where $ can be decomposed as $(ii, f,Z) = 
f{u, v) + g{Z), for some functions /, g and Z stands for mn non-negative variables Zij, i = 1, . . . , m, 
j = 1, . . . ,n. Moreover, Z satisfies some prescribed marginal distributions, i.e., 

E,Z^j=a^(n), EiZij=Pj{v), Zij>0, (14) 

for some functions ai{u) and I3j{v), which satisfy ^^ aj(u) = ^ • f3j{v) = 1. In this section, we will 
only consider u,v such that ai{u) ^ ^ /^jiv) for all i,j, since the remaining cases correspond to 
boundary conditions in our setup and these are treated more easily in the specific instance of the 
maximization (see the relevant Lemmata 12 and 15). Finally, the function gCZ) is an entropy-like 
function and has the form 

ff (Z) = Ei Ej {Zij In Mij - Zij In Z,,) , (15) 

where M S ^"^x" jg ^ matrix with non-negative entries. 

We will use another interpretation of the function ^(Z), which will be handy when we perform 
integrations: let Z denote the variables Z^j for i = 1, . . . ,m — 1, j = 1, ... ,n — 1. Using the 
equations in (14), we may write fl'(Z) as a function of the variables Z . Note that the region (14) 
must be adjusted to account for the non-negativity of the substituted variables. Other than this, 
the "new" function has exactly the same behaviour as ^(Z) but in a full-dimensional space. Thus 
we will also think oi g{Z) as a function of the variables Z*' where the Zij with i = m or j = n are 
just shorthands for the expressions of these variables when substituted by the equations in (14). We 
will refer to this reformulation as the full dimensional representation of fl'(Z). We define similarly 
the full dimensional representation of ^{u,v, Z). 

It is immediate to see that (7(Z) is strictly concave in the convex region (14), hence it attains its 
maximum over the region at a unique point. A tedious but straightforward optimization of (/(Z) 
using Lagrange multipliers, yields: 

Lemma 9. Fix u, v such that ai{u) ^ ^ /^ji^) for all i,j. Let Z* = argmaxz fi'(Z) (the maximum 
taken over the region (14)j. Then 

Zlj = MijRiCj, (16) 

where Ri and Cj are positive real numbers, which satisfy 

Ri ^ Mij Cj = ai and Cj ^ Jvhj Ri = Pj. (1 7) 

3 i 

Moreover, the full dimensional representation of g(Z) decays quadratically in a sufficiently small 
ball around Z*. 

The proof of Lemma 9 is given in Section 9.2. It is relevant to point out that while the values 
of Zj*- are unique in the region (14), the values of Ri,Cj which satisfy (17) are not unique since 
a positive solution {{Ri}i,{Cj}j} of (17) gives rise to a positive solution {{tRi}i,{Cj/t}j} for 
arbitrary t > 0. Modulo this mapping, the uniqueness of Z?- easily implies that the positive Ri, Cj 
which satisfy (17) are otherwise unique. 

Plugging (16) in g(Z), we obtain 

max5(Z) = -E. Y.M,R,Cj HRiC,) := g*{u, v). (18) 



In our setup, it will be hard in general to solve for the Ri, Cj, so that we need to study -^ 
using implicit differentiation. The following lemma will be very useful. 
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dg* dg* 
du ' dv 



Lemma 10. For g* = maxz5'(Z), we have 

du ^ ^ du "^ "^ dv ^—^ ^^'" dv 



Ei..w)^»^|^=-EHc,;^ 



In particular, 

^il±^ = ^-yin(RA^ and ^il±l^ = ^ - Y ln(C ■)^ 
du du ^ ' du dv dv ^ ^ dv ' 

where the Ri,Cj are positive solutions of (17). 
The proof of Lemma 10 is given in Section 9.2. 

5.2 The Logarithm of the First Moment - Maximum 

In this section, we prove Lemma 1. As observed in Sections 4.3 and 5.1, it suffices to study the 
maxima of (f)i{a,P). Lemma 1 is immediately implied by the following lemmas. 

Lemma 11. The only critical points of (j)i{a, f3) are: 

1. (a,/3) = (p^,p^) and (a,/3) = {p*,p*), when Bi,B2,X lie in the non-uniqueness region. 

2. (a, /3) = {p*,p*), when Bi,B2,X lie in the uniqueness region. 

Lemma 12. The function (j)i{a,f3) does not have any local maximum on the boundary of the region 

{(a,/3)|0< a < 1, 0< /3< l}. 
Lemma 13. Let X* = argmaxx *l>i(p''',p~,X) and X° = argmaxx <^i(p*,p*,X). 

1. When Bi,B2,X lie in the non-uniqueness region, the function $i(a,/3,X) has a local maxi- 
mum at the point (p'^ ,p~ ,'X.*) and a saddle point at {p* ,p* ,'K.°). 

2. When Bi,B2,X lie in the uniqueness region, the function <I>i(a,/3,X) has a local maximum 
at{p*,p*,X°). 

We defer the proofs of Lemmata 12 and 13 to Section 9.3. We give the more interesting proof 
of Lemma 11. 

Proof of Lemma 11. We apply Lemma 10 with 

Z -^ X, u -^ a, w -^ /3, ai ^ a, 02 -^ I — a /3i -/^ /3, /32 -^ 1 — /3, 
/^(Q + /3)lnA + (A-l)/i, g^Agi. 

The matrix M is given by 

Lemma 10 gives that the critical points of maxx 'l'i(a,/3,X) must satisfy 
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where 



Ri{BiCi + C2) = a, 
R2{Ci + B2C2) = l-a, 



CiiBiRi + R2) = (3, 

C2{Rl+B2R2) = l-P. 



(21) 



For r = R1/R2 and c = C1/C2 we have 



Bic+l 

\ + B2 

Bir + l 



a 



1 — a 

/3 



a+i/(A-i) 



(l/A) 



1/{A-1) 



r + B2 1-/3 



J+1/(A-1) 



(l/A) 



1/(A-1) 



Bic+1 
C + B2 

Bir + 1 

r + B2 



A-l 



A-1 



The right hand side equations are exactly the equations (3) and in hght of Lemma 7 it must be the 



case that in the non-uniqueness region either that r 



TTq^^ 



or r = c = -^^ , . By (6), we 



■ l-gT - ■ - i-q- 

hence obtain that the critical points for the first moment in the non-uniqueness region are given 
by (a,/3) = {p^,p^) or (q,/3) = {p*,p*). Similarly for the uniqueness region. D 

5.3 The Logarithm of the Second Moment - Maximum 

In this section, we prove Lemma 2. As established in Sections 4.3 and 5.1, it suffices to study the 
maxima of (?!)2(7, S). Lemma 2 is immediately implied by the following lemmas. 

Lemma 14. Under the hypotheses of Theorem 3, for (a,/3) = (p ,p^), the only critical points of 
(j)2{'y,S) satisfy 7 = 0^, 6 = (3'^. 

Lemma 15. For (a,/3) = (p , p^), the function <p2{'y,6) does not have any local maximum on the 
boundary of the region 

{(7,<5)| 7> 0, a-7> 0, 1 - 2a + 7 > 0, 5>0, f3-S>0, l-2f3 + 5>0}. 

Lemma 16. For (a,/3) = {p^,p~), let Y* = argmaxv '&2(a^)/3^j Y). The function <^2(7,<^,Y) 
has a local maximum at (a^,/3^,Y*). 

The proofs of Lemmata 15 and 16 are valid in general, so that to prove Lemmata 3 and 4, one 
needs to obtain the analogs of Lemma 14. The proofs of these analogs are given in Sections 9.5 
and 9.6, see also Section 5.4. We defer the proofs of Lemmata 15 and 16 to Section 9.4 and give 
here the more interesting proof of Lemma 14. 

Proof of Lemma 14- We apply once more Lemma 10 with 

Z -^ Y, u -^ 7, v -^ (5, ai -^ 7, 02 = 03 -^ a — 7, 04 •(— 1 — 2q + 7 
/3i ^ 7, /32 = /33 ^ /3 - 5, /34 ^ 1 - 2/3 + 5, / ^ (A - l)/2, g ^ Ag2. 

The matrix M is given by 

/ Bf Bi Bi I \ 

Bi B1B2 1 B2 

Bi 1 B1B2 B2 

V 1 B2 B2 Bl J 

We obtain that the critical points of maxY ^2{l^ <^) Y) must satisfy 



M 



dj 



In 



Xl-2a + 7) 
. (a -7)2 



A-l 



\R1R4) 



0. 



(22) 
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86 



In 



5{1 - 2(3 + 6)\^~^ fC2C3'^^' 



(/3 - 6Y 



C1C4 



0. 



where the Ri,Cj satisfy the following equations (symbol o denotes the Hadamard product). 



(23) 



/ Ri\ 




/ c\\ 




/ 7 \ 




/Ci \ 




( Ri\ 




( ^ \ 


R2 


M 


C2 




a — 7 




C2 


M 


R2 




P-5 


Rz 




C3 




a — 7 




C'3 




R3 




(3-6 


\R, J 




\cj 




\ 1 - 2a + 7 / 




\cj 




\rJ 




\l-2p + 5 j 



(24) 



We will write out an explicit form of the equations after establishing the following claim. 
Claim 17. It holds that R2 = R3 and C2 = C3. 

Proof of Claim 17. Observe that the matrix M remains invariant upon interchanging its 2nd, 3rd 
rows and 2nd, 3rd columns, i.e., 

(i?i,i?2,i?3,^4,Ci,C2,C3,C4) satisfy (24) ^^ (iii, i23, ^^2, ^^4, Ci,C3, C2, C4) satisfy (24). 

Since 52 (Y) is strictly concave, this yields the claim. D 

The rest of the proof will work as follows: as in the first moment, we first do some manipulations 
with the equations so that we arrive to a nice form. This form coincides with an inequality version 
of the tree equations (3) for the ferromagnetic Ising model without external field, which can be 
analyzed easily. Let us begin by using Claim 17, to write the equalities (24) as 



Ri{BlCi + 2BiC2 + Ci) 



Ci{BfRi+2BiR2 + R4) 



S, 



R2{BiCi + {B1B2 + 1)C2 + ^2^4) = a - 7, C2{BiRi + {B1B2 + l)i?2 + B2Ri) =(3-5, (25) 



R^Ci + 2B2C2 + BiC\) = 1 - 2a + 7, 



C^iRi + 2B2R2 + B^R^) = 1 - 2/3 + J. 



d"f 



85 



In 



In 



7(1 -2a + 7) 
(a -7)2 

5{l-2p + 5) 
{(3-5y 



A-l 



r2 X AN 

JXn 



\ R\ R4 

A-l / ^2 \ A'^ 



0, 



CIC4 



(26) 
(27) 



We now set ri = Ri/R2,r4 = R4/R2, c\ = C1/C2, C4 = C4/C2. After dividing the appropriate pairs 
of equations in (25), we obtain 



n 



c\ 



7 



^ici + (B^B2 + 1) + B2C4 
a--i Bfci + 2Bi + C4 

5 Bm + {B1B2 + 1) + B2n 



(3-5 Bin + 2Bi + n 

Equations (26) and (27) become 



ta 



C4 



1 - 2a + 7 Bici + {B1B2 + 1) + B2C4 

a - 7 ci + 2^2 + Bid 

1-2^ + 5 Bm + {B1B2 + 1) + B2r4 



(3-5 



{nr^y 



7(1 -2a + 7) 
(a -7)2 



A-l 



(C1C4) 



n + 2^2 + Blr4 



5{l -2^ + 5) 
{p-5f 



A-l 



(28) 
(29) 



(30) 
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Observe that 

{Bin + 2Bi + uXn + 2B2 + B^u) = {Bm + {B1B2 + 1) + B2uf + (1 - BiBifinu - 1). 

For convenience, we now work with d = A — 1. Multiplying the equations in (28), we obtain 
^ 7(1 -2a + 7) jBiCi + (^1^2 + 1) + B2Cif 

'"'"'^ (a -7)2 ■(SiCi + (Bii?2 + l)+S2C4)2 + (l-i?li32)2(ciC4-l)- ^ ^ 

It will be convenient at this point to work with d = A — 1. We plug the first equation in (30) into 
(31), yielding 

^^^' {Bin + (51^2 + 1) + i?2C4)2 ' ^"^ {Bm + {B1B2 + 1) + i?2r4)2 ' ^ ^ 

The second equality in (32) follows by a completely symmetric argument for ci,C4. Observe that 
(32) implies that one of the following three cases can hold. 

Case I : rir^ = c\c^ = 1, Case II : rir^ > 1, C1C4 > 1, Case III : rir^ < 1, C1C4 < 1. 

Case I reduces to 7 = a^, 6 = /3'^ hy (30). Thus we may focus on Cases II and III. We further 
restrict our attention to Case II, with Case III being completely analogous. Use the AM-GM 
inequalities Bici + -B2C4 > 2^/BiB2CiC4, Biri + i?2'"4 > 2\JBiB2rir^ to obtain 

(nr4)^/-^-l< i^-BrB2?i^.c,-l) (,^,^)i/. _ ^ < i^ - B.B2?{nr.- I) 

{BiB2 + l + 2VBiB2CiCi) {BiB2 + l + 2y/BiB2ririY 

It is straightforward to verify the identity 

(51^2 + 1 + 2z/b^) + (1 - BiB2f{z'^ - 1) = ((1 + BiB2)z + 2/b^^ 

Set X = y/rirl, y = yjc\c^. Using the identity with z = x,y and taking square roots, we obtain 

^i/d < (1 + BiB2)y + 2/Bl^ ya ^ {l + BiB2)x + 2./B^2 .33. 

" 2,/WB^y + B1B2 + 1 ' ^ " 2y/WB^x + B1B2 + I ' 



Using the substitution B' = {1 + B1B2) /2\/ B1B2, this can be rewritten in the form 

,,V^<^. ,V.<^. (34) 

y + B' X + B' 

The astute reader will immediately realize the analogy of (34) with the tree equations (3) for the 
Ising model without external field. Moreover, since i?' > 1 we are in the case of the ferromagnetic 
Ising model. For the ferromagnetic Ising model on the infinite tree T^+i, we have uniqueness when 
1 < B < -^^- Hence intuitively the lemma holds when B' < ■^^. However, in (34) we have an 
inequality version of the tree recursions (3) and therefore a bit more work is required. 

In our setting, we obtain that (34) cannot hold when x > l,y > 1. To see this, multiply the 
inequalities in (34) to obtain 

.VV/'^<^^.^^. (35) 

^ - x + B' y + B' ^ ' 

Note that (35) is reversed in CASE III. The following lemma implies that (35) cannot hold in 
Case II, and similarly for the reverse inequality (35) in Case III. 
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Lemma 18. Let d > 2. When B' < l/BdTd+i) = ^, for z > 1 it holds that z^/'^ > ^ "t ^ . 
The inequality is reversed for z < 1. 

We give the proof of Lemma 18 after observing that B' < ^^ reduces to \JBiB2 > r^ ■ This 
concludes the proof of Lemma 14. D 

Proof of Lemma 18. The case z < 1 reduces to the case z > 1 using the inversion z ^ 1/z. Thus 
we focus on proving the Lemma for z > 1. Let f{w) := w'^^^ — B'w'^ + B'w — 1. The lemma reduces 
to proving f{w) > when w > 1 (under the substitution w = z^' ). We have 

f'{w) = {d + l)w'^ - dB'w'^-^ + B'. (36) 

f"{w) = dw'^-^ {{d + l)w -{d- 1)B') . (37) 

The assumptions imply that {d + 1) — {d — 1)B' > 0. Hence, when w > 1, we have f"{w) > /"(I). 
Note that /"(I) > only if {d+l) - {d- 1)B' > 0. Thus, f"{w) is strictly positive when zt; > 1 
and hence f'{w) is strictly increasing. It follows that for tf > 1, we have f'{w) > /'(I) and 

f'{l) = {d + l)-id-l)B'>0. 

Thus, / is strictly increasing for w > 1, yielding f{w) > /(I) = 0. D 

5.4 Some Remarks 

Here we give the idea behind the arguments establishing Lemmata 3, 4. As observed in Section 5.3, 
these lemmata boil down to proving that, for (a,/3) = {p'^,p~), the only critical point of (j)2il,S) 
is at (7,(5) = (a^,/?^) under the respective hypotheses. 

Our point of departure is once again the equations in (32), albeit we proceed with a more 

thorough analysis. Looking at the argument for the case rT] < VB1B2, the point of the proof 
which is subject to tighter analysis is the use of the AM-GM inequalities. These are weak if the 
ratios ri/r^, C1/C4 are much bigger than B2/B1. In Section 9.5, we turn this weakness into our favor. 
Namely, for the antiferromagnetic Ising model without external field, in the region < i? < V^~ , 
the densities p^ are heavily biased towards 1 and respectively, and this reflects at the values of 
ri,r4^,ci,C4^. This observation can be turned into an elementary proof of Lemma 3. For the proof 
of Lemma 4, it is easier to get an explicit handle on the values p^ (and thus on the bias of the 
ri,Cj), which can be used to analyze the equations (32) (see Section 9.6). 

Further, let us make some small minor observations. The reader may have noticed that the proof 
of Lemma 14 applies in a more general setup than the one stated in Condition 1. Namely, in the 
region \/BiB2 > /-T] 1 it verifies Condition 1 for all values of a, /? and not just for (a, (3) = {p^ , p^) • 
The proof also extends in the uniqueness region of the antiferromagnetic Ising model (even with 
external field). While we do not draw this path here, since we are only interested in the values of 
(a, /3) which maximize Eg[Z^ ], one can readily determine explicitly a broader range of a, /3 for 
which Lemma 6 continues to hold. 

Finally, for the reader who may feel that the reduction of Lemma 14 to the tree recursions 
is somewhat artificial, let us clear these doubts: in the early stages of this work we analyzed the 
second moment equations for the antiferromagnetic Ising model without external field at the critical 
point, similarly to [21]. As pointed out in the end of Section 3, one then needs to analyze the second 
moment for a = /3 = 1/2. For these special values of a, /3, the solution of the system of equations in 
(25) can then be penetrated analytically and a much more complicated argument leads eventually 
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to an equality version of the inequality system (34). It should be clear that the method presented 
in the proof of Lemma 14 not only leads to stronger results, but also is an approach that one could 
hope to apply in even more complex settings. 

6 NP-Hardness Results 

In this section, we give the proofs for our NP-hardness results. We will prove the following Theorem, 
which allows us to prove Theorems 1, 2, 3. 

Theorem 5. For a given A, assume that {Bi,B2,\) lie in the non-uniqueness region. Moreover, 
assume that Condition 1 holds. Then, unless NP= RP, there is no FPRAS for the partition function 
of the 2-spin model with parameters Bi,B2,X in graphs with maximum degree A. 

With Theorem 5 at hand, it is straightforward to prove Theorems 1, 2, 3. 

Proofs of Theorems 1, 2, 3. Theorem 3 follows immediately from Lemma 2 and Theorem 5. Theo- 
rem 1 was known to hold for A = 3 and A > 6 ( [24] , [7] ) . Lemma 4 in combination with Theorem 5 
prove the cases A = 4, 5 as well. 

We now prove Theorem 2. In the case of Ising model with no external field, note that Lemma 2 
and Theorem 5 give hardness for A > 3 when 

l^zl <B < ^^—^ where d = A-l. (38) 

Vd + 1~ d + 1 

The theorem will be proved once we show that this can be extended to the region < B < rT , ■ 
Observe that hardness for B, A gives hardness for i?, A + 1. It is easy to check that for d > 2, it 
holds that 

VdTT-1 d- 1 d 

V^TT+1 ^ c^+l ^ d + 2' 
so that the intervals corresponding to A and A + 1 are overlapping for A > 3. Thus it suffices to 
check NP hardness for d = 2 and < i3 < r-^ . This follows from Lemma 3 and Theorem 5, thus 

Vd+l 

completing the proof. D 

We thus focus on proving Theorem 5. Sly's reduction [24], while presented for the particular 
case of the hard-core model, is general enough to cover general 2-spin models. The only part in 
Sly's reduction which needs modifications is in establishing the properties of the gadget he uses. 
Let us first describe the gadget used in [24]. 

Assume m' = 0{n^i'^). The construction has two parts: first, a random bipartite graph G with 
n + m' vertices on each side is constructed. Let the vertices on the i-th side be Wi U C/j, where 
\Wi\ = n and \Ui\ = m' . The edges of the graph are the union of A — 1 random perfect matchings 
between Wi U Ui and W2 U U2, together with a perfect matching between Wi and W2. Thus in the 
resulting graph vertices in W := Wi U W2 have degree A and the vertices in U := Ui L) U2 have 
degree A — 1. We denote the graph distribution defined by this construction as Q. 

The second part of the construction appends k disjoint (A — l)-ary trees of even depth i (with 
£ = 0(log n)) appropriately to G. Specifically, the leaves of each tree form a subset of one of Ui, U2 
and the leaves of the trees form a partition of the set U1UU2 (assuming that |f7i| = \U2\ = k{A — 1) 
for an integer k). The trees have otherwise no other common vertices with the graph G. Denote 
by Ri (resp. R2) the roots of those trees whose leaves are a subset of Ui (resp. U2). 
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Denote by H the final graph. Note that the cardinahties |C^i|, |C/^2| should be chosen appropri- 
ately so that the construction of the graph H is possible. Moreover, without going into the exact 
details, k and i should be chosen with care so that the reduction goes through. This is done in [24] 
and it is enough for us to note that the selection of parameters there can be extended to the 2-spin 
model as well. We should emphasize that the crucial property we need for extending Sly's scheme 
is that the sets Ut have much lower cardinality than the sets Wj. 

Analogously to [24], let P denote the phase of a configuration on H, i.e. P = 1 if the number 
of vertices assigned — 1 in Wi is greater than the number of vertices assigned — 1 in W2 , otherwise 
P = 2. The property of the gadget H which is required by the reduction in [24] is the following: 
in the Gibbs distribution of the general 2 spin model, conditioned on the phase P, the spins of 
the vertices in R := Ri L) R2 should behave roughly as a properly defined product measure on 
{ — 1,+!}^. Moreover, the marginal probability of a vertex in Ri having the spin — should be 
different than the respective quantity of a vertex in i?2- At an intuitive level, this property allows 
H to function as a boolean gadget (depending on whether P = 1 or P = 2). 

We now sketch how the required property of H is established in [24]. This is done in two steps: 

1. First prove that in the graph G, conditioned on the phase P, the spins of the vertices in 
Ui U U2 behave in expectation as a properly defined product measure on {—1,+!}^^'^^^. 
This should not come as a surprise: we have already seen in section 2 that for a A-regular 
random bipartite graph, with high probability over the choice of the graph, in the Gibbs 
distribution, the total weight of configurations from the set T,^ '^ dominate exponentially 
over the rest. An analogous phenomenon occurs for a graph G ^ Q, albeit with a slightly 
different quantitative structure. 

2. Translate the results to R. This is accomplished in a crafty way in [24]: the product measures 
on Ui U U2 turn out that feel like even and odd boundaries in the infinite (A — l)-ary tree, 
providing for an elegant translation of the results. Moreover, the extremality of the measures 
implies that the roots of the trees, conditioned on the phase, are strongly concentrated around 
their expected value. The latter part follows after studying non-reconstruction in the extremal 
measures and proving that the correlation of the spins in R to the spins in U decays doubly 
exponentially in i. 

In Sections 6.1 and 6.2, we give the Lemmas in [24] which require adjustment in the 2-spin model. 
Their proofs follow the same way as in [24] , by comparing the quantities in the graph distribution Q 
with those in the distribution Q. Finally, in Section 6.3, we give a sketch of the proof of Theorem 5, 
which is basically checking that Sly's reduction goes through. 

6.1 The Partition Function for Graphs Sampled from Q 

In this section, we give an outline of the Lemmas which are needed to extend Sly's reduction for 
general 2-spin models. 

Let r] be an assignment of spins to U, i.e. 77 G { — 1; +1} • For i = 1,2, denote by rjf the number 
of vertices in Ui which are assigned spin +1 and let r/^ = m' — -qf (the number of vertices in Ui that 
are assigned spin —1). For a graph G ^ Q, let Z-Q{rj) be the partition function over configurations 
which agree with r/ on U, i.e. ZQ{r]) is the total weight of configurations in cr G {— 1,+!} '•'^•' 
subject to a\u = rj. We have the following analogue of [24, Lemma 3.1]. 

Lemma 19. For any a, /3 and all r] € {0, 1} we have 



(1 + 0(1)) ^^^ — ' ^ ^\ ^'_ -Bl' "^^ A^i +^^2 



Eg[Z^'''] \{a- x*)"!! {I- a- (3 + x*)'^i-'^2 {/3 - x*)'^2 
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where x* is the (unique) non-negative solution of the equation BiB2{a — x){f3 — x) = x{l — a — f3+x). 
In particular, when (a, /?) = {p^,p^), it holds that 

^ G ^^'^ = (l + o(l)) - ('o+)''i"fl - a'^)'>Ua~)''^ (1 -(7")''2" 

Eg[Z^'^] ^'^"^'>\l-q+r(l-q-)m'^1 ) [i q ) [q ) [i q ) , 

where 

\B2 + Q+ + Q- + BiQ+Q-) ■ ^ ' 

The proof of Lemma 19 is given in Section 7.2.1. Seeking high probabihty lower bounds on Z^ [t]), 
we study the second moment of Z^ [j]). 

Lemma 20 (analogue of Lemma 3.5 in [24]). When Condition 1 is valid, for (a, /3) = {p^,p~) and 
all rj G {0, l}*^ we have 



Mi^G^f] ^ ^^^ ^ y yB2 + Q~J j \ \B2 + 

where Q~,Q'^ are given by (2) and C* is given by (39). 

The proof of Lemma 20 is given in Section 7.2.2. Once again, it turns out that the second moment 
method fails to give high probability bounds, as is captured by the following lemma. 

Lemma 21. When Condition 1 is valid, for (a,/5) = {p'^,p~) it holds that 

hm ^g[(^g"^^^)'] = (1 _ ^2)-(A-l)/2 .^ _ i)2^2)-l/2 

7°'- 

'g 
where uj is given by (4). 



M^li^v)] 



The proof of Lemma 21 is given in Section 7.2.3. We then apply the small graph conditioning 
method to obtain the following analogue of [24, Lemma 3.9]. 

Lemma 22. For A > 3, suppose that for Bi,B2,X in the non-uniqueness region. Condition 1 is 
true. Then, for {a, (3) = {p'^,p~) asymptotically almost surely over the choice of the graph G ^ Q, 
for every r] £ {—1,1}'^ , it holds that 

The proof of Lemma 22 is given in Section 8.2. 

6.2 The Reconstruction Problem in the Extremal Measures 

The final ingredient we need is proving that a strong form of non-reconstruction holds in the 
extremal measures. Let us first introduce some notation. 

Denote by Tp^n^/s. the subtree of Ta rooted at p up to depth L Let Sp^i be the leaves of Tp^n^^. 
Denote by ^p,^,+ the marginal probability that the root p is occupied in a configuration a generated 
by the following process: 
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1. Sample a configuration a from the projection of the measure /i"^ on Tp^i^j\, 

2. Conditioning on tlie configuration as on Sp^i, sample a configuration a from the Gibbs dis- 
tribution on Tp^i^/\. 

Note that the configuration of the vertices in Sp^i is a random vector, so that ^p,^,+ is a random 
variable. Define similarly Xp^i_. We need the following analogue of [24, Lemma 4.2]. 

Lemma 23. Assume that A > 3 and Bi,B2,X lie in the non-uniqueness region. There exist 
constants Ci(^5 Bi,B2, X), C2(^7 Bi,B2,X) > 0, such that for all sufficiently large i, 

^[\Xp,i,± -q^\> exp(-Ci^)] < exp(-exp(C2^)). 

We omit the proof of Lemma 23: Martinelli et al. [17, Proof of Theorem 5.1] prove the required 
strong concentration for the variables -^p,£,± for the Ising model with arbitrary external field (in a 
different form than the one stated here but completely equivalent). Their results can be translated 
to the general 2-spin model, after carrying out the reduction described in [23]. Finally, in the case 
of hard constraints Lemma 23 follows by the respective result in [24, Lemma 4.2] (slightly extended 
in [7, Lemma 30]). Of course, one could also apply the methods in [17] or [24] to directly get the 
result in general. 

6.3 Proof of Theorem 5 

The reduction of [24, Section 2.2] goes through as is, once we establish [24, Theorem 2.1]. Thus, 
we focus on checking the relevant parts of its proof. 

Lemma 1 together with Lemma 19 suffice to prove [24, Lemma 3.3]. Lemma 1 together with 
Lemma 22 suffice to prove [24, Theorem 3.10]. The proof of [24, Lemma 4.3] is clearly valid for the 
2-spin model. Having these ingredients, the proof of [24, Theorem 2.1] pulls through in the 2-spin 
model as well. 

7 Calculating Moments 

7.1 Random Bipartite A- Regular Graphs 

We use the information we have obtained for the first and second moments to prove Lemma 5. Its 
proof is a very technical calculation. 

Recall that Q^ = , ^ j, . To aid the presentation, it will be convenient to also have the following 
expressions defined beforehand. 

El := BiQ-Q+ + Q-+Q+ + B2, E2 := B1Q-Q+ + BiB2{Q+ + Q") + B2, 
E3 := (1 - BiB2fQ^Q+ + {l + Bi{Q+ + Q-) + BfQ-Q+^ (^B^ + B2{Q~ + Q+) + Q~Q+) . 

(41) 
These expressions will emerge in the analysis of the asymptotics of the first and second moments, 
as is captured by the following lemmas. We are also going to need the explicit values of X* = 
arg maxx 51 (X) and Y* = argmaxx52(Y) when a = p'^,l3 = p~,7 = oP',5 = /3^. These are given 
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by: 






I vli y*i2 yh y*ii \ 

^ ^ ^ ^ 

2/21 ^22 2/23 2/24 

2/31 2^32 2/33 2/34 

2/41 2^42 2^43 2/44 / 



1 



/ Bf{Q-)HQ+r B,Q-{Q+r B,Q-{Q+f {Q+f \ 

Bi{Q-fQ+ B1B2Q-Q+ Q-Q+ B2Q+ 

Bi{Q-fQ+ Q^Q+ BiB2Q~Q+ B2Q+ 

\ {Q-? B2Q- B2Q~ Bl ) 



(42) 



(43) 



The easiest way to argue about the vahdity of these values is to use the strict concavity of the 
functions gi (X) , (^2 (Y) and just check that the required equahties for the critical points are satisfied. 

Lemma 24. For {a,P) = {p^,p^), we have that 



n^oo ie^^iC^'/S.x*) 27r^ ' \ E? 

Lemma 25. For (a,/3) = {p~^,p~), we have that 
,. Eg[(Zg-^)^] 1 m^2(A-i) 



{Q+Q-)^ElE,Y''^' ({Q+Q-fE2E^\^'^ (l - (A - 1)V)^'/' 



where 7* = a^ ,5* = /3^, Y* zs given by (43) anc/ a; ^y (4)- 

The proofs of Lemmas 24 and 25 are given in Sections 7.1.1 and 7.1.2 respectively. We need to 
state one more lemma before we can give the proof of Lemma 5. 

Lemma 26. It holds that 2$i(q,/3,X*) = $2(a^/3^Y*) 

Proof of Lemma 26. It clearly suffices to prove 

/2(a^/32) ^ 2/1 (a, /3) and 52(Y*) = 2gi(X*). 

Both of these can be verified by straightforward brute force algebra. A neater way to do this, using 
more directly the uncorrelated property that Condition 1 represents, is to use entropy arguments. 
We illustrate this for convenience in the case of /i and /2, with the same technique extending to 
51,52 as weU. 

Define A, B to be independent Bernoulli variables with success probabilities a, j3 respectively. 
It is straightforward to see that /i(a,/3) = —(^H{A) + H{B)^ where H{X) is the entropy of the 
random variable X. Let A' = (^1, ^2)1 where Ai,A2 are independent copies of the random variable 
A. Define similarly B'. Using the definition of entropy, it is easy to check that f2{oP,P'^) = 
-{H{A') + H{B')). But clearly H{A') = H{Ai) + H{A2) = 2H{A) and similarly for B,B', thus 
verifying that /2(a^/3^) = 2/i(a,/3). 

The observation which extends this idea for 52 ( Y* ) = 2gi (X* ) is that the solution Y* , as given 
in (43), is the tensor product of X*, as given in (42), with itself. We omit the tedious details. 
We should emphasize that the proof of Lemma 26 does not rely on the fact that Q^,Q~ are a 
function of the critical probabilities on the infinite tree, i.e., we do not need Lemma 7, but rather 
the representation given in (6) and the uncorrelated property that Condition 1 represents. D 
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Proof of Lemma 5. Employ Lemmata 24, 25, 26 to obtain 



Eg 
L := lim 



G 



(Zen 



Eg[Z^'^- 



(1 - (A - 1)V)-^/^ (a/3(l-a)(l-/3))^-^ ^(Q^Q^^EEsy/^^A/, 



„ ,(Q+Q-)4i?|i?3^ fQ+Q-E2Y' {Q+Q-fE2E^ 

K 



'G 
where 

u-fiOl^ , 

El^ J \ El J El 

Thus 



L = (1 - (A - l)V)-^/^(a/3(l - a)(l - /3))^"i (* 



:i-(A-i)v -^^ 



2 2^-l/2r 1 {Q+Q-fE^E,---^''-'^" 



a2/32(l_a)2(l_^)2 ^7 

By (6), we have a = ^^(^+fQ') and P = ^lil^^, so that 



(i?2 + Q-)HB2 + Q+)2(l + BiQ-)2(l + SiQ+)2 

The lemma follows after observing the easy to prove identity 

E1E2E3 2 

I — to . 



{B2 + Q")2(S2 + Q+)2(l + i?ig-)2(l + i?iQ+)2 

D 

7.1.1 The asymptotics of the first moment 

In this section we prove Lemma 24, i.e. compute the asymptotics of Eg[Zg ] for (a,/3) = {p~^,p^), 
see the relevant Sections 4.1 and 4.3. 

Proof of Lemm,a 24- We will use the full dimensional representation (see Section 5.1) of gi(X.): 
while we are still going to refer to the values x*,, in what follows the only free variable is going to 
be xii. By the approximation 



(^^\ = (1 I 0(1)) ^ / ^ ^n(blnb-alna-(a-b)ln(a~b)) 

\anj V2^V a{b-a) 



we obtain 



-1/2 1 



^ = (l + o(l))--(a/3(l-a)(l-/3))' 



lgn$i(a,/3,X*) ^ ^ ^^27r 



E(n-«)d^ -"<"■'-•-'■'""" 
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Note that in the above sum terms with |j;ii — a;|^| < 0(1) dominate exponentiahy over the rest and 
in a small enough ball around x^^ the function decays quadratically, as established by Lemma 9. 
Thus, standard integration arguments, displayed in full rigor in [13, Chapter 9], yield 

n^oo ie^^iC"'/''^ ) 27rV ^ ^^ V V-Li ^-'Z 

exp Tr-2-(2;ii) ^ "a;ii 



Let -ff := ^-^(2;^]^). Since 51 is strictly concave and it achieves its maximum at xli we have H < 0. 
It follows that 



i/:-(a«''^^^^^)--'-'-"' 



and consequently 



„S"»i|iSS^ - ^M' -»)(!- «)"^-'"X(-«'n4^"^'^ 



Xi 



We have 

+ + -7; + 



dxf^ V^ii ' OL — x\\ ' fi — x\\ ' 1 — a — /3 + xii 

Substituting the values of a,/3 by (6) and the solution for x\y by (42), we obtain 

S2 + Q- + Q+ + BiQ-Q+){B2 + SiS2Q~ + 51^2(3+ + BiQ-Q- 



-H 



n 



B1B2Q-Q+ 
BiB2{Q+Q-f 



so that 



^5 J 






D 



7.1.2 The Asymptotics of the Second Moment 

In this section we prove Lemma 25, i.e., compute the asymptotics of Eg[(Z^ )^] for {a,j3) = 
{p~^iP~)- We have 

7,5 

where Eg[l^' ] was determined in Section 4.2 (see also the relevant Section 4.3). 

Proof of Lemma 25. As in the proof of Lemma 24, in the following we are going to refer to the full 
dimensional representations (see Section 5.1) of 32, ^2, i-e. our free variables are going to be {yij} 
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with 1 < i, j < 3, so that the rest of the yij are just shorthand for their respective values (in terms 
of the free variables) . We have that 

n^ 7,(5 ^^ ^ 

The function $2(7, (5, Y) has a unique maximum at (7*,(5*,Y*). Thus, in the above sum, terms 
with II (7, 5, Y) — (7* , 6* , Y*) II < 0(1) dominate exponentially over the rest. Moreover, in the region 
II (7, 6, Y) — (7* , 6* , Y*) II < 0(1) the function <1>2 decays quadratically. As in the proof of Lemma 24, 
we thus obtain 

li- i^'1S7L = (^) '' (7*^*(l - 2a + 7*)(1 " 2/3 + ^*)(« - 7*)^(/3 - S*f^ ^ 



A 



n 

A 
4 \ 2 



00 fOO 
ij 



1 / \2(A-1) , __ 

^(a/3(l-a)(l-/3)) H y,^ 



A 
4 \ 2 



i)7:/:(7;^/:-'"'"'''"'"^^ ^^* <-> 



where H denotes the Hessian matrix of <1>2 evaluated at (7*, (5*, Y*) scaled by 1/A and the operator • 
stands for matrix multiplication. 

We thus focus on computing the integral in (44). This will be achieved by two successive 
Gaussian integrations. First, we split the exponent of the quantity inside the integral as follows. 

3 3 

i(7,<5, Y) • H • {j,5,Yf = 1(7,5) • H^,^ • (7, 5)^ - ^Y • (-Hy) • Y^ + ^E^^^^^i' 

i=i j=i 

where H^^^ denotes the principal minor of H corresponding to 7, 5, Hy denotes the principal minor 
of H corresponding to Y and 

1 / d^<^2 , , d^^. 



^•' A V ' djdyij ' " d6dy. 
It follows that 



-00 / \J—oo 

Let T denote the row vector with entries Tjj. Note that Hy is the Hessian of (72 (Y) evaluated 
at Y*. Since 52 (Y) is concave and §2 achieves its maximum at Y*, we have that Hy is negative 
definite. Hence, 



A 

fOO 



_J /""^ gi(7,5,Y).H.(7,5,Y)T^Y 1 ^ ^#(t.(-Hy)-^-T^+(7.^)-H^.^-(7.^)^ 

V2^f J-00 J (Det(-HY))^/' 



24 



We are thus left with the task of computing the integral 



)^ fOO fOO ^ , 

/ / e#(T-(-HY)-i-T^+(7,5)-H,,,-(7,<5r)_ (45) 

J — 00 J— 00 

It can easily be seen that the expression in the exponent is a quadratic polynomial in 7, 6, so that 
we may again perform Gaussian integration. We have 

^ (T • (-Hy)-^ • T^ + (7, <5) • U^,s ■ (7, Sf) = Dj' + Ejd + F5\ (46) 

for some complicated D,E,F, which can be determined explicitly by (46), once we perform the 
substitutions. Hence, 

provided that ADF — E'^ > 0, so that the integration is meaningful. We will check this condition 
later. 

Combining equations (44), (45), (47), we obtain 

_ A 

1 / \2(A-1) 

= ^(a/3(l-a)(l-/3)j ■ R, (48) 

where 

R := (Det(-HY) fl y*A [ADF - E^r^'\ 

At this point, we resort to Maple to perform the necessary computations. The solution in the form 
of (43) is particularly handy, since it rationalizes the expressions in the computations. We have: 

Note that by Lemma 8, we have ADF — E^ > 0, so that the integration in (47) is indeed meaningful. 
After plugging the above values into R and substituting back in (48), the lemma follows. D 

7.2 Gadget 

In this subsection we prove Lemmas 19 and 20. We will need the following bound. 

Lemma 27 (general form of Lemma 3.2 in [24]). Let bi,...,bi and yi,...,y£ be non-negative 
integers. Let a = bi + ■ ■ ■ + b£ and x = yi + ■ ■ ■ + y^. Assume yf < b-i for each i G {1, ...,£}. Then 
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Proof. We have 



The lemma now follows from the fact that x'^/a < Vi/hi + • • • + yf/b'j (which follows from the 
Cauchy-Schwarz inequality). D 

7.2.1 First Moment 

Proof of Lemma 19. By (7), recall that 



E,[.-]^A.-^.«(^;)(;)^|:^^^^.s- 



where 



^a,P,x 



fan\ f(i-—a)n\ 

\xn) \{l3-x)n' j^x 0(^—01— P+x)n 



Similarly, we have 



BfB, 



A-l 



\x<rain{a+ri~ /n,l3+-q~ /n} / \x<min{a,/3} / 



where 



'' MR-x^n+r,::) 



n+m' 



G ^ " ( n+m' \ ^ 

\Rn+r,~) 



From Lemma 27 and the fact that rjf,?]-^ ,r]2,r]2 < n^'^, we have the following estimate on the 
ratio of K^^'^{r]) and k^'^. 

G '" =(l + o(l)) ^^ — ' ^ ^\ '^ B^i "2^ (49) 

It is easy to see that for any a, /3, the logarithm of the function Kq'^ scaled by n is exactly 
the full dimensional representation of the function 51 (X) (see Section 4.3). Since fl'i(X) is strictly 
concave with respect to X, it has a unique maximum. With these observations, it can easily be 
calculated that the value of x = xu for which (71 (X) achieves its maximum is given by the unique 
positive solution of 

BiB2{a - x*)(/3 -x*) = x*{l-a-l3 + x*). 
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By Lemma 9, the logarithm of Kq'^ decays quadratically in the distance from x* . Let 

A = {x: \x-x*\ < n"^/^}. 

The relative contribution of terms outside A to Eg[Z^ ] is exp(— ri(n^'^)) and hence can be 
omitted. Similarly, using (49), the relative contribution of terms outside A to Eg[Z^ (ry)] is 

exp(— fi(n^'^)) and hence can also be omitted. 
For x&Awe have 

^''^(^) ,i\M)) a^^ (1 - art pv^l - p)4 ^,._,-^ 

k'^^''' (a - x*)*?! (1 - a - /3 + x*rt-V2 (/3 - x*)^2 ^ 

Thus, 



Eg[Z°''^] \{a-x*r^{l-a-P + x*rt-'^HP-x*r^ ^ ' 



A-1 



which proves the first part of the lemma. 

When (a, /3) = {p'^,p~), we express a, (3 and x* = x^i in terms of Q^,Q^ using equations (6) 
and (42). We obtain 

a^i (1 - art (3^^ (1 - Pr^ ^r,t~4 _ 



(a - x*r^ (1 - Q - /3 + x*)^i -'^2 (/3 - x*r^ 



B. 



2 



^2(l-«)(l-/3) \'"' f a{l-a-(3 + x*) y^ /^a^^^^^x;^\ ^^ 
l-a-;5 + x* y V52(a-x*)(l-a)y VS2(/3 - x*)(l - /3), 



(52 + <3+)(52 + Q~) \'" /l + BiQ-\"^ fl + BiQ 



''h /I I R, n+ \ ''2 



,-B2 + Q+ + Q- + 5iQ+Q-y V^2 + Q-y V^2 + (3^ 

and hence 

^[^5^<">i - a .„„)c- (. {^4^rr u i'-p^rr . m 



where C* is given by (39). Using the fact that Q^,Q^ satisfy (3), we have 



B2 + Q- J " \B2 + Q+ 

Hence, we obtain 



Eg[Z,"'^l 



(l + o(l))C*(l + Q+)"*(l + Q-)^ 



G J 

Q+ N'^i' / 1 yt / Q- ^2" / 1 ^2+ 



l + Q+J \l + Q+J \l + Q-J \l + Q 

The second part of the lemma follows after observing that q^ = tttw j ^ — Q^ = , , q± ■ D 
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7.2.2 Second Moment 

Proof of Lemma 20. By (10), we have 



A 



where 

,,Y^A/ L,n \'f R,n \ / n \ 

,~,i Vyii^' y^2?^, ?/i3f^, mW -^^^ Vyii^i y2jra, ys^-n, y4jW \V\\n, yun, ..., yiinj 
y, T^{'2yii+yi2+yi3+y2i+y22+y3i+ys3)nT^{y22+y24+y33+y34+y42+y43+'2y44)n 

Li=7, L2 = L3 = Q-7, L4 = l-2a + 7, 
i?^ = 5, R^ = R^ = j3-5^ i?4 = 1 - 2/3 + (5, 

and Y stands for the nonnegative variables y^j, 1 < «,J < 4 such that 

Ej-yij = Li and ^liyij = ^j- 
Similarly, we have that 



where 



LiTi \ -i-r / -Ro-n \ 1 1 n + m 



^ ,_i Kyu'iT'^ yi2n, Visn, ViAnJ f^ \yijn, y2jn, y^jU, yijiij \ynn, y^n, ..., y^iU 

y 0(2^11+2/12+^13+2/21+^22+^31+2/33)™ o(y22+3/24+J/33+j/34+3/42+y43+2j/44)n 

-£i=7 + ^r' 1/2 = -£3 = 0-7; -£4 = 1 -2a + 7 + r?j'", 
iii = (5 + 7?^, R2 = R3 = /3-5, R4 = 1-2/3 + 6 + 1]+, 

and Y stands for the nonnegative variables yij, 1 < i,j < 4 such that 

Ej-yij = i'i and Xlj-yy = Rj- 

For a given Y, consider Y such that yij = ijij for all i,j G {1,2,3,4}, except 2/14 = yu + ?/f, 
2/41 = ?/4i + f?^ , and ^44 = 2/44 + 771' — ?/j~ — r/^ . This is always possible, except at the boundary, but 
as we will see shortly, these cases can be safely ignored. By Lemma 27, we have 

!^1^ = (1 + od))^"- <^ - ^° + -<r-'-" (^-v+ ")"' gw-.,.-.,.-^ 

'^G yllyilvAA 



Let 



^ = {(7,<5, Y) : ||(7,<5, Y) - (7*,<5*, Y*)||2 < n~'/'}, 
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where -y* = a'^, S* = f3'^ and Y* is given by (43). By Condition 1 and Lemma 9, the relative contri- 
bution of terms outside A to Eg;[(Z^ ) ] is exp(— r2(n^/^)) and hence can be omitted. Similarly, 
using (53), the relative contribution of terms outside A to Eg[(Z^ (ry)) ] is ex.p{—Q{n^''^)) and 
hence can also be omitted. For (7, J, Y) G ^ we have 

^t^ (yJjV(yl^)^.-(ylJ-'-^r-'y.- 2 

Using equations (6) and (43), we express a,/3 and Y in terms of Q~^,Q^. For (7, 5, Y) £ A, we 
obtain that 

^''""(ri) _^^^ ^ ^^^^^ fBlil - 2a + 7*)(1 - 2/3 + 5* 



Kg V y44 



Vi / x*„,* \ V2 



Yvu X" f s*yl4 



BlvUl -2a + r)J \Blyl,{l -2(5 + 6* 

^ ^^''\B2 + Q+ + Q- + BiQ+Q-) \B2 + Q- ) KB^ + Q+J ' 
Plugging the last equation into (51) and (52) we obtain the lemma. D 

7.2.3 Proof of Lemma 21 

Proof of Lemma 21. By equations (40) and (50), we have 



(%[^5''W])' "-°°(lEg[^G'1)'' 



G 
The lemma follows, after using Lemma 5. D 

8 Asymptotically Almost Surely results 

8.1 Overview 

In this section, we apply the small graph conditioning method to obtain Lemmata 6 and 22. The 
following theorem is taken from [21]. The notation [X]m refers to the ?7i-th order falling factorial 
of the variable X. 

Theorem 6. Let Aj > and 5i > —1 be real numbers for i = 1, 2, . . .. Let r(n) — )■ and suppose 
that for each n there are random variables Xi = Xi{n), i = 1, 2, . . . and Y = Y{n), all defined on 
the same probability space Q = Qn such that Xi is nonnegative integer valued, Y is nonnegative and 
E[y] > (for n sufficiently large). Suppose furthermore that 

1. For each k > 1, the variables Xi,. . . ^X^ are asymptotically independent Poisson random 
variables with EXj — )■ \i, 

2. for every finite sequence r?ii, . . . , m^ of nonnegative integers, 

"^''^''''^E[y]"™"^^ ^n(^^(i + ^^))"'- (54) 
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3. E.M'<oo. 

I E[y2]/(E[y])2 < exp ( Y.i ^i^i) + o(l) asn^oo. 

Then Y > r(n)E[y] asymptotically almost surely. 

Recall that Zq is the measure of configurations in S^ for a random A-regular bipartite graph 
G. Let Xi be the number of circles of length i in G. Clearly Xj = iff i is odd. We have the 
following lemmas. 

Lemma 28 (Lemma 7.3 in [21]). Condition 1 of Theorem 6 holds for even i with 

_ r{A,i) 

■^i — : ) 

I 

where r{A,i) is the number of ways one can properly edge color a cycle of length i with A colors. 
The proof of Lemma 28 is given in [21]. 
Lemma 29. For (a,/3) = {p ,p^), we have that 

^Q[^G ^^\ ^x,{l + 5i) and ^^ ^ V^ ^ ^A,(l + ^,), 
^g[Z-/] %[^f'^(r/)] 

where 

(l-5iS2)^/2(Q+Q-)V2 



(I + SiQ+)(l + BiQ-){B2 + Q+){B2 + Q- 
The proof of Lemma 29 is given in Section 8.2. 
Lemma 30. It holds that 



i/2' 



exp Y. AAM=(1--^)"^^"'^/'(1-(A-1)V 



^even «>2 

where uj is given by (4). 

The proof of Lemma 30 is given in Section 8.2. 

Lemma 31. For (a, /3) = {p^,p'^) and for every finite sequence mi, . . . , ni]^ of nonnegative integers, 
it holds that 

% [Zq {r])[X2]rrn • • • [-^2fc]mfc] , TT /' \ n ^ X W^' 

a,/3 >ll[Mi- + di)) , 

^Gi^G \ i=l 



and 

^- "^^ V'/;K>-2Jmi • • • L-'^2feJmfc 



Eg[zf{ v)[X2U--- [X2k' ' " 



11(^^(1 + ^^))"- 



The proof of Lemma 31 is identical to [21, Proof of Lemma 7.5]. 
We are now in position to prove Lemmata 6 and 22. 

Proofs of Lemmata 6 and 22. We verify the conditions of Theorem 6 in each case: Condition 1 
holds (Lemma 28), Condition 2 holds (Lemma 31), Conditions 3 and 4 hold (Lemma 30 and 
Lemmata 5, 21). Lemmata 6 and 22 follow by picking r{n) = 1/n and r{n) = l/\/n respectively. D 
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8.2 Random Bipartite A-regular Graphs 

In this section, we prove Lemma 29. We give the proof for the variables Zq and omit the proof 



for the variables Z^ {?]), which can be carried out the same way as in [24, Lemma 3.8] and along 



the lines of Section 7.2.1. 

Recall that r(A, i) denotes the number of proper A-edge colorings of a cycle of length i. 



Lemma 32. r{A,i) = (A - 1)^ + (-l)^(A - 1). 

The (simple) proof of Lemma 32 is given in [21, Lemma 7.6]. 

Proof of Lemma 29. For sets 5, T such that S C Vi, T C V2, |5| = an, \T\ = j3n, denote by Ys^t the 
measure of the configuration a that S, T induce, i.e. for a vertex v E V{G), a{v) = 1 iff u G 5 U T. 
Fix a specific pair of S, T. By symmetry, 

nZc'^X,] ^ E[Ys,tX,] 

Eiz^'"] nYs,T] ■ ^ ' 

We now decompose Xi as follows: 

• ^ will denote a proper A-edge colored, rooted and oriented i-cycle (r(A,i) possibilities), in 
which the vertices are colored with red(i?), blue(i?) and white(VK) such that no two adjacent 
vertices in the cycle have both R or B color, i.e. the only assignments which are prohibited 
for an edge are (R, R), (B, B). The color of the edges will prescribe which of the A perfect 
matchings an edge of a (potential) cycle will belong to. R denotes vertices belonging to 5, B 
denotes vertices belonging to T and W the remaining vertices. 

• Given ^, C, denotes a position that an z-cycle can be in (i.e. the exact vertices it traverses, in 
order) such that prescription of the vertex colors of ^ is satisfied. 

• 1^^(^ is the indicator function whether a cycle specified by ^, C, is present in the graph G. 

Note that each possible cycle corresponds to exactly 2i different configurations ^ (the number of 
ways to root and orient the cycle). For each of those ^, the respective sets of configurations C, are 
the same. Hence, we may write 



2i 



Let Pi := Pr[l5,^ = 1]. It follows that 



? c 

In light of (55), we need to study the ratio E[1s'^t|1^^(^ = 1]/E[15^t]- At this point, to simplify 
notation, we may assume that ^, C are fixed. 
We have shown in Section 4.1 that 

•sr-^ \Sn) \{/3-S)n) ^fa^(l-a-/3+^)n 

where 5 denotes the number of edges between 5", T in one matching. 
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To calculate E[l5^T|lg^(^ = 1], we need some notation. For colors ci, C2 G {R, B, W}, we say that 
an edge is of type {ci,C2} if its endpoints have colors ci,C2. Let ji,J2 denote the number of red 
and blue vertices in the coloring prescribed by ^. For A; = 1, . . . , A, let ai{k) denote the number of 
edges of color k of type {R, B}, 02 (A;) denote the number of edges of color k of type {W, W}, di{k) 
denote the number of edges of color k of type {R, W}, d2{k) denote the number of edges of color k 
of type {B, W}. Finally, for j = 1, 2 let aj = ^^ aj{k) and dj = Y^^ dj{k). By considering the sum 
of the degrees of R vertices, the sum of the degrees of B vertices and the total number of edges of 
the cycle, we obtain the following equalities. 

ai + di = 2ji, ai + d2 = 2J2, ai + 02 + di + ^2 = i- (56) 

A straightforward modification of the arguments in Section 4.1 yields 



^ / fan~ai{k)—di{k)\ / {l—a)n—a2{k)—d2{k) \ 

I \(l3-Si.)n-a-i{k)-d2(k)) 

)] 

fc=l \ 5fc V I3n-ai{k)-d2{k) J 



¥,]Y ll = 1] = \(°+/3)" TT i \^ '^'='" ^^il^-^k)n-ai{k)-d2{k)J 5kn+ai(k) rj{l-a-l3+5k)n+ai{k) 

[ S,T\ £,,C — \ — J_J_ I / ^ ,n~ai{k)-a2{k)~di{k)-d2{k)\ 1 2 



By Lemma 27, we have 

(an-ai{k)-di{k)\ , _ \ ai(fc)+di(A:) 



I OLn\ 
Wn) 

( {i--a)n-a2{k)-d2(,k) \ . a2{k)+d2ik) / ox \ ai(k)+d2{k) 

y{0-Sk)n-ai{k)-d2(k)) (l-a-P + 6k\ ^' ^' P - Ok \^' ^' 



( (^-a)n \ \ I- a J \l- a- 13 + 6, 



3fc 



/n—ai{k)—a2{k)-di(k)—d2{k)\ 

\ I3n-ai(k)-d2{k) I ^ oai{k)+d2{k) ^-^ _ n\a'^(k)+di{k) 

\l3n) 

By the same line of arguments as in the proofs of Lemmata 19, 20, we obtain that 
¥.[Ys^t\M,c = 1] {5*Y^~''^{a - 5*)"i+'^i(/3 - 5*Y^^'^^{BiB2Y' 

^\Ys,t\ '" Q;'^i+di^ai+d2(l _ Q,)a2+d2(l _ /3)a2+di ' 

Clearly Pi ~ ?i~* and for given ^, the number of possible C, is asymptotic to 

aJ'i/3J2(l - q)*/2-J'i(1 - /3)^/2-i2n\ 
Thus, for the given ^, we have 

X;^PiE[Fs,t|1^,C = 1] (1 - a - /3 + 5*Y^-''^{a - 5*Y^+'^^{I3 - 5*Y^+'^KBiB2Y^ 

^\Yst\ aai+di-ji j^ai+d2-J2(\ _ Q,)a2+a!2+il-i/2n _ ^^a2+c^l+j2-«/2 

Utilize (56) to express everything in terms of i,ji,J2,ci2 to get 

(l-a-/3 + 5*)*-2ji-2i2(a - 5*)2ii(/3 _ S*Y^^BiB2Y^~'^^^'^^^^ 
ah(3J2{l - a)i-Ji(l - P)^-^^ 

Substituting 

Q+{l + BiQ-) ^ Q~{l + BiQ+) ^, BiQ~Q 



R. 



B2 + Q- + Q+ + B1Q-Q+" B2 + Q-+Q+ + B1Q-Q+' B2 + Q~+Q+ + BiQ-Q+' 
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it follows that 



R 



Bl^ -i+2ii+2J2 5«2 (g+ )il (Q- ) i2 



We thus obtain 



where 



(1 + BiQ-)h{i + BiQ+)n{B2 + Q-y/'^-n{B2 + Q+y/^-i^ ' 
E[Ys,t\H( = 1] 



X 



V 



Ws,t\ 

BfQ+{B2 + Q 



x'y^^z^^w"^, 



Thus, by (55), we obtain 

E[YXi] r{A,i] 



1 + BiQ- 



1 + BiQ+ 



-, w = B1B2. 



E{Y] 



2i 



/ J ^jl,J2,a2y ^ ^ J 



(57) 



JlJ2,a2 



where ajjj Jj^qj is the number of possible ^ with ji red vertices, j2 blue vertices and 02 edges of type 
{VF, ly}. To analyze such sums, an elegant technique (given in [12]) is to define an appropriate 
transition matrix. The powers of the matrix count the (multiplicative) weight of walks in the 
implicit graph. In our setting, the transition matrix is given by 



A 



f y y \ 

z z 

1 w 

\ 1 w J 



Note that a rooted oriented cycle with specification ^, ^ corresponds to a unique closed walk in the 
implicit graph defined by A. We have that Tr(^') = X^fci^i) where ej are the eigenvalues of A 
(j = 1, . . . , 4). We thus obtain 



E[Y] 



2i 



j=i 



A computation gives 

ei = -62 = Bi(l- B1B2 



Q+Q- 



^ {1 + BiQ+){l + BiQ-) 
It follows that for even i, we have 



, 63 = -64 = Bi^{B2 + Q+)iB2 + Q-). 



EjYXj] _r(A,i) 
E[Y] i 



1 + 



{l-BiB2)y^Q+Q'~ 



V(l + BiQ+){l + BiQ-){B2 + Q+){B2 + Q- 

thus proving the lemma. 

Proof of Lemma 30. Using Lemma 32, we have 



D 



E A.<^^= E 



r(A,i) 



• UJ 



E 



(A - 1)* + (A - 1) 



• OJ 



even i>2 



even i>2 



even i>2 
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Observe that XIiM %T" ~ "^l^l-*- ~ ^^) ^^^ ^^^ 1^1 < 1- Clearly this holds for u and by Lemma 8, 
it also holds for (A — 1)l(J. It follows that 



^ AA' = -i(ln(l-(A-l)V) + (A-l)ln(l 



a;2 



even i>2 
thus proving the lemma. D 

9 Remaining Proofs 

9.1 Proof of Lemma 8 

In this section, we give the proof of Lemma 8. We prove only the inequality for w, since the 
inequality for uj* is the condition for non-uniqueness, see the discussion in Section 3. We will in 
fact establish a slightly stronger inequality, as is captured by the following lemma. To aid the 
presentation, let us work with x = Q^,y = Q~ , d = A — 1. 



Lemma 33. For any d> 2, it holds that Bixy + BiB2{x + y) + B2 > {d — 1){1 — BiB2)^/xy. 

Proof of Lemma 8. Let W := Bixy + BiB2{x + y) + i?2- Observe 

W = {Bix + l)(y + B2) - (1 - BiB2)y = (B^y + l)(x + B2) - (1 - BiB2)x. 

Using the two expressions of W and the AM-GM inequality, we obtain 

W^ = {Bix + l)(y + B2){Biy + l){x + B2) + (1 - BiB2fxy 

- (1 - BiB2){y{Bix + l)iB2 + y) + x{Biy + 1)(52 + x)) 

< {ViBix + l)(5iy + l)iB2 + x)iB2 + y) - (1 - ^iSs)^^)'. (58) 

Trivially 

{Bix + l){Biy + 1){B2 + x){B2 + y) > xy > {I - BiB2fxy, 

so (58) and Lemma 33 give 



v/(Bix + l)(y + B2)iBiy + l)(x + B2) - (1 - ^1^2)^^ > (d - 1)(1 - B1B2). 

which after trivial manipulations reduces to the desired inequality. D 

Proof of Lemma 33. The x, y satisfy 

jB^y + lY , .fB^x + 1 

x = A \ — — and y = A 



y + B2 J \x + B2 

It follows that 

^(?T^) =^(7!^) ^^{iB^^ + ^)iy + B2)Y = y[{Bw + l){x + B2)] . (59) 
Let W := Bixy + B1B2X + BiB2y + -62- Observe that 

{Bix + l){y + B2) = Bixy + B1B2X + y + B2 = W + {l- BiB2)y, 
{Biy + l)(x + B2) = Bixy + BiB2y + x + B2 = W + {l- BiB2)x. 
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Hence, (59) gives 

x{W + {l- BiB2)yt = y{W+{l- B,B2)xf . 



Expanding using Newton's formula gives 

d / ,N d 

X 

k=0 "^"^ fc=0 

which is equivalent to 



I — n V / u—r^ V / 



fc=2 

Since x ^ y^ this can be rewritten as 



W\x -y)=xyJ2 (f) ^'"'(1 - BiB 

I — o \ ■' 



2nx'~'-y'~' 



W' = xyj: ['^W'-Hl - B,B,f ( ""^:f" ) . (60) 



^k-l_ 


.yk-l 


X - 


-y 


x'^-1- 


_yk^l 



Claim 34. For k>2 andx,y >0 with x j^y, it holds that > (k - l)(xy)(''~^)/^. 

x-y 

The simple proof of Claim 34 is given at the end. Using Claim 34, (60) gives 

d 



W'>Y1 (f) (^ - ^W^~\'^ - BiB2f{xyY 



\k/2 
k=2 ^ ^ 

or equivalently 

W^ + Y, (^) W'^-'^il - BiB2f{xyfl^ > J2 (^) kW^'^il - BiB2f{xyfl^ . 

k=2 ^ ^ fe=2 ^ ^ 

^ v ' ^ v ' 

C D 

Using again Newton's formula and the identity (^) = ^(/.Zi)) we have 

C={W + {1- B^B2),/^f - dW'-^l - i?iB2)Viy, 

D = d{l- B,B2)V^Y it I 1) ^'"'(^ - BiB2)'~HV^)''-' 



(61) 



fc=2 

= d{l - BiB2)V^[{W + (1 - BiB2)V^)''~^ - W-' 

Thus, (61) gives 

W>{d-l){l-BiB2)V^, 

which is exactly the inequality we wanted. Finally, we give the proof of Claim 34. 
Proof of Claim 34- Since k >2, observe that 

"^'"^ ~ y'~' = x^-2 + x'^-^y + ...+ /-2 >{k- l)('(xy)('=-^)('=-2)/2y^^'"'^ = {k- l)(xy)('=-2)/2. 
X — y V / 

The inequality is an application of the AM-GM inequality to x^~'^, . . . ,y^~'^. Equality holds iff 
X = y. D 

This completes the proof of Lemma 33. D 
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9.2 Proofs for Section 5.1 

Proof of Lemma 9. We first give the proof when Mjj > for all i,j. Then it is easy to see that the 
function 5'(Z) attains its maximum in the interior of the region (14), since 

InMj,' — In Zij — 1 —7- +00 as Z^ — t- . 



... ..^jj iiL^lj ± —r -rv^ cxo ^jj 



Since (7(Z) is differentiable in the interior of the region (14) (by ctj 7^ 7^ /3j), we may formulate 
the maximization of g'(Z) using Lagrange multipliers. Let 

A(Z, A, V) = g{T) + Ai(V Z,^ - a^) + A' ( V Z,,- - /3,). 



i 



The optimal Z must satisfy 
5A 



dZ,, 



I'd 



^^ In Mij - In Zij - 1 + Aj + A^- ^^ Z^j = exp (Aj + A^- - l) Mij . 



Setting Ri = exp(Aj),Cj = exp(A'- — 1) one obtains that the stationary points of g{Z) have the 
form stated in Lemma 9. The equalities (17) come from the conditions ^ = M^ = 0. 

We next treat the case where some entries of Mij are zero. Recall the convention (Section 4.3) 
that InO = —00 and OlnO = 0. To see why these conventions are relevant, notice that whenever 
Mij = 0, g(Z) can have a maximum at Z* only if Z*- = 0. Going back to the interpretations of 
5f(Z) in the first and second moments, the reader will immediately notice that these are in complete 
accordance with the combinatorial meaning of the variables Zij. Let P = {{i,j) ■ Mij = 0}. We 
may restrict our attention for the maximum of g(Z) in the region 



J2j Zij = ai{u), for i G 1, . . . , m, 

(62) 



Y^i^ij = Pj{v), for j G l,...,?i 



Zij>0, for(i,i)^P, 

Zij=0 for(i,i)GP. 

The proof given above (under the assumptions M/j > for all i,j) now goes through for the region 
(62), as it can readily be checked. Since Mij = for {i,j) G P, the form of Z* stated in Lemma 9 
is consistent with the region (62). 

Finally, the desired quadratic decay is an immediate consequence of the stationarity of Z* , the 
strict concavity of g and Taylor's expansion. D 

Proof of Lemma 10. From (17) we have 

-i?.)j;M,,C,+i?.5]M,,-C, = -a., (63) 

-CAY,M,,R, + C,Y.M,,-R. = 0. (64) 



Then 

4 = 1 j = l \ ^ 

d 



-R^ (C, \n{Ri) + C, \n{Cj)) + (^C^ {R^ HR^) + Ri HCj))+ 



jrRi I Cj + Ri 
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which, using (63) and (64), simphfies to 

dg* ^, ,^,1 ( d \^^^ ^ r.l\^.. d 



Y^HR^ (-rA Y,m,,c, + r E^^^a^^. 

fl^<^^^ [{k^) XJ^%^. + C, (^EM,,|-i?,,^^ 

j = l \ ^ ^ j = l 4=1 / 





n 



d 

1=1 



u 



9.3 Proofs for Section 5.2 



Proof of Lemma 12. To avoid unnecessary coniphcations, we give the proof for B1B2 > 0. The 
case with Bi = (essentiaUy hard-core model) is covered in [6, Claim 2.2]. By (19) and (20) we 
have 

where the Ri,Cj satisfy (21). By (21), we have 

a _Ri B1C1 + C2 (i _Ci BiRi + R2 

l-a~R^'Ci + B2C2 ' 1-P~C^'Ri + B2R2 ■ 

da ' df3 



Hence, the derivatives — — , — — - become 



da \ a \Ci + B2C2 

9^^, / 1-/3 f BiRi + R2 Y 
dp ""V ^ \Ri + B2R2) 

The following claim, which we prove at the end, will be helpful. 

Claim 35. When Bi > 0, -B2 > 0, the function 

Bix + y 



' infA^f^i^^^rv m 



(68) 



f{x,y) ■■-- 



X + B2y ' 



defined for x,y > with xy ^ is upper and lower bounded by strictly positive constants depending 
only on Bi, B2- 
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Observe that Claim 35 is applicable to our case since not all of Ri , R2 or Ci , C2 can be equal 
to (by equations (21)). 

With Claim 35, we can easily analyze the boundary cases. We consider the boundaries with 
a(l — a) — )• 0"'" and prove that the derivatives with respect to a go the right way. The remaining 
cases (and their combinations) are completely analogous. We consider cases. 

(a) a —^ 0~^ . Then 1 — a — )■ 1. By Claim 35 and (67), we obtain -^ — )• +00. 

(b) 1 — a —7- 0"^. Then a — t- 1. By Claim 35 and (67), we obtain -^ — t- —00. 

We now give the proof of Claim 35. 

Proof of Claim 35. Observe that for t > it holds that f{x, y) = f{tx, ty), so that we may assume 
x + y = 1. Thus, Bix + y and x + B2y are convex combinations of {1, i?i} and {1, B2} respectively. 
Since we assumed that i?i,i?2 > 0, f{x,y) is bounded by strictly positive constants. D 

This concludes the proof of Lemma 12. D 

Proof of Lemma 13. Denote by Hi the Hessian of ^i evaluated at {p'^ ,p~ ,X*). We keep only 
a,(3,xii as free variables in $1 (the rest of the Xij are substituted by (8)). We check that the 
principal minors along the diagonal are positive, the first principal minor being the lower right 
corner of the matrix. This is equivalent to applying Sylvester's criterion in a rearranged version of 
—Hi. We use Maple to perform the necessary computations and plug into the resulting formulas 
the values of p , X* as given in (6) and (42) respectively. 

Denote by Pi the determinant of the i-th leading principal minor. We obtain 

AEiE2 _ AEKB2 + Q~){l + BiQ-){l + iA-l)u;) _ AEfjl - {A - l)^u) 

' B,B2Q-Q+' ' BiB2{Q-rQ+ ' ' BiB2{Q+Q-f ' 

where Ei,E2 and uj are given by (41) and (4). Trivially Pi,P2 > 0. The fact that P3 > is just 
Lemma 8. 

To prove that (p* ,p* ,X°) is a saddle point, observe that the respective expressions for Pi, P2, P3 
can be obtained by substituting Q~ ^ Q*,Q'^ ^ Q* , where Q* = jr^. We have that P3 < by 
Lemma 8. D 

9.4 Proofs for Section 5.3 

Proof of Lemma 15. To avoid unnecessary complications, we give the proof for B1B2 > 0. The 
case with Pi = (essentially hard-core model) is covered in [7, Lemma 12]. 
By equations (26) and (27), we have 

7(l-2a + 7) \^"^ / Pj Y) = 1,, ( («-7)^ / 7(l-2a + 7) ^j Y] 



d-f \\ (a -7)2 J \R1R4J / l7(l-2a + 7)V (a - 7)^ PiP4y 



^2 




5{l-2f3 + 5)\^^^ f Cl \\_. { {^-^? /^5(l-2/3 + (5) Cl 



86 y\ (/3-(5)2 J \C1C4J J y6{l-2p + 5)\ {^ - 5f Cid 

where the Ri,Cj satisfy (25). By (25), we have 

7(1 - 2a + 7) _ P1P4 {BfCi + 2P1C2 + C4)(Ci + 2P2C2 + BJCa) 

(a -7)2 ~ Pi ■ (PiC7i + (PiP2 + l)C2 + P2C4)2 
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6{1 -2/3 + 5) _ C1C4 (BfRi + 2B1R2 + Ri){Ri + 2B2R2 + B^R^) 



{13 - 5)2 d {BiRi + (S1S2 + l)R2 + B2Rif 



Hence, the derivatives — — , -— ^ become 

(77 00 



In 



(a - 7)2 / {BlCi + 2B1C2 + Ci){Ci + 2B2C2 + Bid] 



dj \^7(1 - 2a + 7) V (^1(^1 + (^1^2 + 1)^2 + ^2^4)2 

^=1 (^ (/g - '^)' A^i'^i + 2^1 ^2 + R4){Ri + 2^2 i?2 + B^Ri) ■ , 
9<5 ''l^<5(l-2/3 + (5) V {BiRi + {BiB2 + l)R2 + B2Ri)^ ' '' ^ ^ 

Before proceeding we need the fohowing claim, which we prove at the end. 
Claim 36. When Bi > 0, i?2 > 0, the function 

{Bfx + 2Biy + z){x + 2B2y + B^z) 



f{x,y,z) :-- 



{B^x + {B^B2 + l)y + B2zY 



defined for x,y,z > with xyz ^ is upper and lower bounded by strictly positive constants 
depending only on Bi, B2. 

Observe that Claim 36 is applicable to our case since not all of Ri, R2,Ra or Ci, C2, C4 can be 
equal to (by equations (25)). 

With Claim 36 at hand, it is immediate to treat the boundary cases. We consider the boundaries 
with 7(0 — 7)(1 — 2q + 7) — )• 0+ and prove that the derivatives with respect to 7 go the right way. 
The other cases (and their combinations) are completely analogous. We consider cases. 

(a) 7 — >■ 0"*". This case may only occur if 1— 2a > 0. Then a— 7 — )■ a > and 1— 2a+7 — )• 1— 2a > 0. 
By Claim 36 and (69), we obtain -^ — ;• +00. 

(b) 1 — 2a + 7 —7- 0*^. This case may only occur if 2a — 1 > 0. Then 7 — )■ 2a — 1 > and 
a — 7— )-l — a>0. By Claim 36 and (69), we obtain -$^ — )■ +00. 

(c) a — 7 — 7- 0+. Then 7 — )• a > and 1 — 2a + 7— t-I — a>0. By Claim 36 and (69), we obtain 

We now give the proof of Claim 36. 

Proof of Claim 36. Observe that for t > it holds that f{x, y, z) = f{tx, ty, tz), so we may assume 
that x + y + z = 1. Thus, Bfx + 2Biy + z, x + 2B2y + B^z and Bix + {B1B2 + l)y-\- B2Z are convex 
combinations of {-Bx, 2Bi, 1}, {1, 2^2, -Bl} and {Bi,BiB2 + 1,-62} respectively. Since we assumed 
that Bi,B2 > 0, f{x,y,z) is bounded by strictly positive constants. D 

This concludes the proof of Lemma 15. D 

Proof of Lemma 16. For (a, j3) = {p'^,p~), denote by H2 the Hessian of $2 evaluated at {a^, (3'^,Y*). 
We keep only ^,6,yij with i,j = 1,2,3 as free variables in $2 (the rest of the yij are substituted 

by (11)). 

We use Sylvester's criterion to check that — H2 is positive definite. We check that the principal 
minors are positive, the first principal minor being the lower right corner of the matrix. This is 
equivalent to applying Sylvester's criterion in a rearranged version of — H2. Note that the first nine 
principal minors are guaranteed to be positive since they correspond to the Hessian of the convex 
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function —g2(Y) evaluated at its minimum point. Thus we need to check only the remaining two 
principal minors. We use Maple to perform the necessary computations and we plug into the 
resulting formulas the values of p^ from (6) and the Y* as given in (43). 

Denote by Pi the determinant of the i-th leading principal minor. We have: 

^^EfEl{B2 + Q~f{l + BiQ-f {l + (A - l)uj^) 

{B,B2Y{Q-f\Q^ 

A9^f E|(l - (A - l)2a;2) 



^10 - /^ „.NQ/^ xl4/^,xl2 



Pii = det(-H2) 



14/^ , \14 



{BrB2f{Q-y\Q+) 



where Ei,E2 and ui are given by (41) and (4). Trivially Pio > 0. And we also have Pn > by 
Lemma 8. D 

9.5 Proof of Lemma 3 

We now prove Condition 1 for the Ising model without external field and for d = 2, that is 
Bi = B2 = B , X = 1 and A = 3. See Section 5.4 for the proof overview. 

Our point of departure is (32). As in the proof of Lemma 2, we may assume that rir4, C1C4 are 
either both greater than 1 or both are less than 1. We are going to argue that (32) cannot hold in 
these cases. For Bi = B2 = B, multiplying the equations in (32) gives 

(Be. + iB^ + 1) + B..nBn + iB^- . 1) + Br.f = " /'C^^ ,' " ^ ^T^'" 7 '' ■ ^ 

{ciCi)^!'^ — 1 (rir4)^/"^ — 1 

We will prove that for < S < \-" , when rir4 7^ 1 (and similarly when C1C4 7^ 1) 

(Bn + (B' + l) + Buf> ^'-''"T^^'-\-'\ (72) 

In the region < -B < \-~ the values of a and /3 are heavily biased towards 1 and respectively, 
and this reflects at the values of ri, r4, ci, C4. The required bias of the values of a, /5 is captured by 
the following lemma. 

Lemma 37. Assume < B < --= . For (a, B) = (p+,p~) it holds that = > -. 

The proof of Lemma 37 is given in Section 9.5.1. We utilize Lemma 37 to prove: 

Lemma 38. Assume < B < y=-— • For (a,/3) = {p~^,p~) andri,r4,ci,C4 satisfying the equations 
in (28), (29), it holds that 

.n>-.^anrfc4>---^. 

4 1 4 1 

• // in addition rir^, C1C4 > 1, then ri > - • -—r and C4 > - • -— r. 

9 is^ 9 i?^ 

The proof of Lemma 38 is given in Section 9.5.2. We now prove inequality (72). Rewrite (72) as 

B^in + r^f + 2B{B^ + l)(ri + r^) + (B^ + 1)^ > (1 - B^f (1 + ,/7^) , 

or 

B\n + r^f + 2B{B^ + l)(ri + r^) + 45^ > (^ _ B^f,/¥^. (73) 

At this point, we split the analysis into the cases rir^ > 1 and rir4 < 1, the first being considerably 
harder. 

40 



Claim 39. Inequality (73) holds when rir4 < 1. 

Proof of Claim 40- We prove the stricter inequality (since 1 > i? > and ri,r4 > 0) 

2Bri > ^/rlr4. 

By item 1 of Lemma 38 and the assumptions, we have 

2 1 , 

25ri > :j • ^ > 1 > Vnn- 



D 



Claim 40. Inequality (73) holds when rir4 > 1. 

Proof of Claim 40. We prove the stricter inequahty (since 1 > B > and ri,r4 > 0) 

B'^{ri+r4f > ^/¥^. (74) 

We first identify the regions where (74) is hard to prove. Note that B'^{rf + 2rir4 + r|) > AB'^rir^, 
so that (74) is true if AB'^rir^ > yjrir^. Thus we may assume that 

> rir^. (75) 



16^4 
Since rir^ > 1, by item 2 of Lemma 38 
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9'B 
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B'rl > j,^, (76) 



so that (74) is true if g^^ > y^rirj. Thus we may assume that 

44 
rm > ^^. (77) 

We are now ready to prove (74). Using (75), (76) and (77) we obtain 

B^irl + 2rir4 + r|) > B^rf + Ir.r^) > ^ [^ + '^ ' ^) > ^ ' \> VnU- 

D 

9.5.1 Proof of Lemma 37 

By the remarks in the end of Section 3, for the Ising model with no external field it holds that 

a _ W 

2 



— - and Q~^Q = 1. For ease of presentation, let x = Q^ , y = Q ■ By (3) and (6), we have 



a x(By + l) , /By + 1\ 

— where x = I — - I and 2; > 1 > y > 0. 



1- a y+B \y+B 

Since xy = 1, we have that 



where -^ = —^ — — and < y < L (78) 



I- a y{y + By y^ \y{y + B 
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Proof of Lemma 31. For < -B < ^^ -*- = 3 — 2-v/2, we want to prove that > - • —-7, which 

V 2+1 1 — a 9 B-^ 

in hght of (78), is equivalent to y < (3/2) ^^S^. For convenience, let c := (3/2)^372. The 

equation for y in (78) may be rewritten as 

(^y + Bf = y{By + if =^ B^y^ + {-I + 2B)y'^ + {I - 2B)y - B^ = Q ^ 

{y - 1) (52y2 + (_b2 j^2B- l)y + B"^) = =^ B'^y^ + (B^ + 2B- l)y + B^ = 0, 

since by (78) we are interested in < y < 1. 

Let f{z) := B'^z'^ + (i?^ + 2B — l)z + B^. By Lemma 3, we know that y is the unique root of 
f{z) for < 2; < 1. Observe that /(O) = -B^ > 0, so the desired inequality will follow if we prove 
that f{B'^c) < 0. Now, / may be rewritten as 

9 / 1 1-2B -B^' 

f{z) = B'z[z + - -^ 

For z = B'^c and < i? < 3 - 2\/2, we have 

1 _ B-^c^ + 1 (3 - 2V2)^c=^ + 1 ^ 
^^~z~ B^c - c 52 ' 

l-2B-B^ 1 - 2(3 -2\/2)- (3-2^/2)2 _ 16\/2- 22 
^2 ^ ^2 - ^2 • 

A straightforward calculation gives ^^ — < f < 16\/2 — 22, thus proving f(B'^c) < 0. The 

lemma follows. D 

9.5.2 Proof of Lemma 38 

We need the following lemma. 

Lemma 41. When B1B2 > and ri,r4,ci,C4 satisfy (28) and (29), it holds that 

a 1 + ri/^2 1 + Bici 13 



1 - a 1 + B2ri ' 1 + C4/51 1 - /3 

We first give the proof of Lemma 38 and then give the proof of Lemma 41. 
Proof of Lemma 38. When Bi = B2 = B, Lemma 41 gives 

a 1 + ri/B l + Bci (3 



1-a 1 + Sr4 ' 1 + C4/5 1 - /3 ■ 

Since r4, ci > 0, it follows that 

where we have used Lemma 37 to bound j^ and -^. This proves the first item of the lemma, 
after observing that | • -g^- — 1 > | • -^ for any B < 1/3. 

To prove the second item of the Lemma, note that when rir^, C1C4 > 1 it holds 

1 + ri/B 1 + ri/B _ri 1 + Bci 1 + B/ca_B 

1 + Sr4 ^ 1 + B/ri ~ B ^" 1 + Ci/B ^ 1 + c^/B ~ 7^' 
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It follows that 

ri> B • and C4 > i? • — -; — . 

1 — a p 

The desired bounds now follow after using Lemma 37 to bound -p^ and — i^. D 

Proof of Lemma 41- It will be convenient to have the following variables at hand. 

Bfci + 2Bi + Ci ci+2B2 + Blci 



BiCi + {BiB2 + l) + B2Ci' BiCi + {BiB2 + l) + B2Ci' 

Bfri + 2Bi+u „ ri + 2B2 + B^ri 

-Ki — — --— — — : — , JriA 



Bm + (51^2 + 1) + ^2^4 ' BiCi + {B1B2 + 1) + B2C4 ■ 

Note that the Ci, C4, i?i, i?4 do not have anything to do with the i?j, Cj in the proof of Lemma 14, 
but this notation is convenient since the equations in (28), (29) may be written as 

nC, = ^-, r4C4 = 1^^^^, (79) 

a — 7 — 7 

5 1-2/3 + 5 



Eliminating 7, 5 from (79) and (80) respectively, we obtain 

a _ 1 + nCi /3 _ 1 + cii?i 

1 - Q ~ 1 + r4C4 ' 1 - /3 ~ 1 + C4i?4 ■ 

We have the following loose bounds on i?i, i?4, Ci, C4. 



(81) 



Bi<Ru ^4<^, Ci<^, B2<Ci. (82) 

i^l 1^2 

These bounds may be proved by picking any of the inequalities, multiplying out and using B1B2 < 1. 
It is immediate now to check that the bounds in Lemma 41 follow after combining (81) and (82). D 

9.6 Proof of Lemma 4 

Lemma 42. For the hard-core model (that is, Bi = 0, B2 = I) for A = 3, 4, 5 and any a, l3 £ (0, 1) 
the solution of (22), (23) and (24) satisfies ^ = o? and 5 = /3^. 

Proof. Let d := A — 1 and let ci := C1/C2, C4 := C4/C2, n := R1/R2, r^ := R4/R2 (the same 
notation as we used in the proof of Lemma 14). Equations in (32), for Bi = and B2 = 1, become 



(ciC4)^/^-l= ;^ and (rir4)^/'^-l= ;;^ (83) 

(l + r4)^ (1 + C4)^ 

Let X = (rir4)-^/'^, y = (ciCi)^''^, a = l/r^, and b = I/C4. We can rewrite equation (83) as follows 



l + 2a + a^x'^ , l + 26 + 62yrf 

W = ;5— and X = ; -^-. (84) 

l + 2a + a2 1 + 26 + 62 ^ ^ 



Solving for a and 6 we obtain 



^^^ and 6 = ^^^-+^, (85) 

x"- — y y°- — X 
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where 

'■2-(y-l)(x'^-l) and ql = {x-l){y''-l). (86) 



ia 



Note that a > and & > 0. Assume y > 1. Then x > 1 (otherwise the sides of (83) would have 
different signs) and also y < x'^ (the right hand side of (84) is a convex combination of 1 and x ) 
and hence \qa\ > y — 1- Thus, using a > in Equation (85), we obtain qa > and similarly q^ > 0. 
Now assume y < 1. Then x < 1 and also y > x'^ (again, the right hand side of (84) is a convex 
combination of 1 and x'^) and hence \qa\ > 1 — y- Thus, using o > in Equation (85), we obtain 
Qa < and similarly ^^ < 0. 

Because of symmetry between x and y we only need to consider two cases (in the case x = y = 1 
we have 7 = 0^ and 6 = 0^): 



• CASE 1: l<y<x<y^,qa:= ^J {y - l){x'^ - I) , Qb := V(x - l)(y'^ - 1), and 

• CASE 2:0<y'^<x<y<l,qa:= -y^iy - l){x<i - 1), q^ := -y^{x - l)iyd - 1). 
From Equation (24) one obtains the following expressions for a and /5 (see equation (81)) 

a _ T^n + 1 ^ P T^ci + 1 



and 



l-« '-^^r, + l 1-/3 ^-^^^c + l 

Solving for a, (3 and using the parametrization with x, y, a, h we have 

y(l + a)2 -l-a + ab 



a 



/? 



x{l + bY + y{l + af + 2ah- l 
x{l + hf -I -h + ab 



(87) 



x(l + 6)2 + y(l + a)2 + 2a6 - r 
Recall that a and /3 satisfy: 

a _ A(l-a-/3)^-i /3 _ A( l - a - /3)^ -^ 

(1-/3)^-1 
and hence 



l-a-/3 


(1- 


q)^-i ' 1- 


Q-/3 






a(l 


-a)^-i=/3(l 


- - /3)^- 


~i 


is equivalent to 


(assuming 07^/?) 








a(l- 


a)^-i-/3(l- 


-/3)^-i 





Equation 

which using the expressions for a and /3 from (87) becomes 

{y{l+af -~l-a+ab){x[l+bf +ab+aY-{x{l+bf -l-b+ab){y(l+af +ab+bY) ^ q /qqn 

The final part of the proof will require some computational assistance (just to manipulate polyno- 
mials). First we are going to plug-in the expression for a and h from (85) into (90). Then we are 
going to use the expressions in (86) to reduce the powers of qa and g^ occurring in the expression 
obtaining an expression of the form 

Coo + ClO^a + CoiQ'b + Ciiqaqb-, 

where coo, coi, cio, cn are polynomials in x and y. Then we are going to show that in CASE 1 we 
have Coo > 0, coi > 0, cio > 0, cn > and in CASE 2 we have coo > 0, coi < 0, cio < 0, cn > (and 
hence (90) cannot be zero) . This will be accomplished by reparameterizing x = {ty + y'^)/{t + 1), 
factoring the expressions, and observing that: 
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• factor y — 1 occurs with even power in cqo and cn, 

• factor y — 1 occurs with odd power in cqi and cio, and 

• all the other factors have all coefficients positive. 

The details of the argument appear in Appendix A. D 
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A Appendix: Computer Assisted Proofs for the Hardcore Model 

A.l Case A = 3 

d = 2; 

F = Factor[((y*(l + a)"2 - 1 - a + a*b)*(x*(l + b)"2 + a*b + a) "d - 

(x*(l + b)'2 - 1 - b + a*b)*(y*(l + a)-2 + a*b + b) "d)] [[3]] ; 

Comment: now F contains the left-hand side of (90) (we took the third factor; the other two 
factors are —1 and the denominator of (90)). 

H = Factor [Expand [F /. 

{a -> (-1 + y + Qa)/(x-d - y) , b -> (-1 + x + Qb)/(y-d - x)}]][[3]]; 

Comment: now H contains the left-hand side of (90) after substituting the values of a,b given 
by (85) (the result has three factors: l/(x — y'^)'^'^~^, l/{x'^ — y)'^'^~^, and the factor we assigned to 
H; note that the first two factors cannot have value zero). 

da = Exponent [H , Qa] ; 
For [i = da, i >= 2, i — , 

H = Expand[H /. {Qa~i -> Qa'd - 2)*(1 - x"d - y + x~d*y)}]]; 

db = Exponent [H , Qb] ; 
For [i = db, i >= 2, i — , 

H = ExpandEH /. {Qb"! -> Qb"(i - 2)*(1 - y'd - x + y~d*x)}]]; 

Comment: now H contains the left-hand side of (90) (multiplied by (x — y^-^''^^^ (^x'^ — y)'^'^^^) 
after reducing the powers of qa and qb using (86). 

cOO = Factor [Coefficient [Coefficient [H, Qa, 0], Qb, 0]]; 

cOl = Factor [Coefficient [Coefficient [H, Qa, 0], Qb, 1]]; 

clO = Factor [Coefficient [Coefficient [H, Qa, 1], Qb, 0]]; 

ell = Factor [Coefficient [Coefficient [H, Qa, 1], Qb, 1]]; 

Comment: the following reparameterization will reveal the signs of cqo, cqi, cio,cii. We show the 
expressions for d = 2; for larger d we will print the first 6 factors and check the positivity of the 
last factor by Mathematica. 

uOO = Factor [Expand [cOO /. {x -> (t*y + y~d )/(t + 1)}]] 

uOl = Factor [Expand [cOl /. {x -> (t*y + y~d )/(t + 1)}]] 

ulO = Factor [Expand [c 10 /. {x -> (t*y + y"d )/(t + 1)}]] 

ull = Factor [Expand [ell /. {x -> (t*y + y"d )/(t + 1)}]] 

OUTPUT: 



47 



- tt^-Tq (-1 + y)^y^(l + t + y)(4:t + 28*2 + 84*3 _^ ^4q^4 _^ -^^q^s _^ 34^6 _^ 28i^ + 4t^ + 12ty 
(1 + t)^ \ 

+ Mt^y + 2A8t^y + iOOt'^y + 380t^y + 212t^y + 6ify + 8t^y + Ay"^ + 50ty^ + 228*^^2 + 532t^y^ 
+ lUt^y"^ + 570t^y2 + 26it^y^ + 64i'^y2 + Gt^y'^ + IGy^ + 141i/ + 502*^^3 + 954t^y^ + 1054t^/ 
+ 683t^y3 + 248**^2/3 + idfy^ + 4*^^^ + 32/ + 233t/ + 694*^/ + lOSQt^/ + 976i^/ + 521t^/ 
+ 179t^/ + 43i^/ + 5t^/ + 40y^ + 243ty^ + 592^2^^ + 759^3^5 _^ gOOi^'y^ + 352t^y^ + 16U^y^ 
+ 43t^2/^ + 4i^y^ + 32y^ + 159ty^ + SlSt'^y^ + 387*^/ + 385tS'^ + 302i^y^ + 139t^y^ + 26t^y^ 
+ t^y^ + IGy'^ + 63ty'^ + 132^2^7 + 241t^y'^ + 322i''y'^ + 236t^y^ + 72^^^"^ + Gfy"^ + 4y^ + 19t?/^ 
+ 82t2y8 + 201*3/ + 234t^y^ + llOt^y^ + Ibt^y^ + Uty^ + 69*^/ + 137*^^9 + lOOt^/ 
+ 20t^?/^ + 10ty^° + 44*^2/10 + 54t3yio + 15t^y^° + 6ty" + IGt^y^^ + Gi^y" + 2ty^2 _^ ^2^12 



' -l + yfy'^{l + t + y)(4t + 24*^ + eot^ + 80t^ + 60t^ + 24i^ + if + Sty + 50t'^y 



{1+tr 

+ 126t^y + 164t^y + 116t^y + 42^*^^ + 6fy + Ay^ + 38ty^ + 136*^^2 + 248*^/ + 254^*^^ _^ 143^5^2 

+ AGt^y^ + 6fy^ + 12y^ + 89ty3 + 258*^^3 _^ 333^3^3 _^ 311^4^3 ^ 133^5^3 _^ 35^6^3 _^ g^7^3 

+ 20/ + 118t/ + 274*2/ + 3i2t3y4 _^ 197^4^4 _^ 92^5^4 _^ 3g^6y4 _^ ^^7^4 _^ 20y5 + g^ty^ + 165*2^5+ 

156t3y5 + 127t^y^ + 87t^y^ + 29t^y^ + 2fy^ + \2y^ + 42t/ + 65*^/ + 94^3^6 ^ 113^4^6 

+ 70t^/ + 12^6/ + 4y^ + Wty'' + 36*2^"^ + 8U^y'^ + 90t''y'^ + 30i^y^ + 6i/ + 35*^/ 

+ 65*3/ _^ 4Q^4y8 ^ g^y9 _^ 25*2/ + 30*3/ ^ 4^yio _^ 12*2^1° + 2ty^^ 



- -y(-l + yfy^{\ + * + y) (4* + 20*2 + 40*3 + 40*^ + 20*^ + 4**^ + 10*y + 50*2y + 96*3y 

(1 + *j v 

+ 88*^y + 38*^y + 6**^^ + 4y2 + 38*y2 + I20*2y2 _^ 174^3^2 _^ i26*''y2 + 44^5^2 _^ g^6y2 _^ ^g^s 

+ 107*/ + 267*2y3 + 324*3/ ^ 197^4^3 ^ 55^5^3 ^ g^6y3 ^ 32^4 ^ ^^^^^4 ^ 371^2^4 ^ 3gg^3y4 

+ 179*^/ + 47*^/ + 5**^/ + 40y^ + 183*y^ + 310*2^^ + 250*3y5 + \\8t^y^ + 28*^y^ + 2t^y^ 

+ 32^*^ + 119*/ + 164*2/ + 132*3y6 + 5gi4y6 ^ g^5^6 ^ ^^yl _^ 47^y7 ^ gg^2y7 ^ ^^t'yl ^ -^2t^y7 

+ 4/ + 13*/ + 23*2/ + 8*3/ + 4*/ + 2*2y9 



- -g (-1 + yTyUt + 20*^ + 40*3 + 40** + 20**^ + 4*'^ + 10*y + 48*^^ + 88*3y + 76**y 

(1 + *j v 

+ 30*^y + Afy + 4y2 + 38*y2 + mt^y^ + 164*3y2 + UQt^y^ + 42*^^2 _^ 7^6^2 _^ ^Qy3 _^ ^^^^yS 

+ 253*2/ + 294*3y3 + 173t4y3 _^ 4g^5y3 _^ g^6y3 _^ 32^4 _^ ^^^^y4 _^ 339^2^4 _^ 310^3^4 _^ 139^4^4 

+ 39*^/ + 4**^/ + 40y^ + 173*y^ + 268*2^^ + 192*3^^ + 92*^y^ + 20*^y^ + 32^*^ + 109*/ + 132*2/ 
+ 102*3/ + 40*^^*^ + 16y^ + 41*^^^ + 54*2^^^ + 40*3y^ + 4y^ + ll*y^ + 20*2^^ + 4*/ 
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A.2 Case A = 4 

d = 3; 

F = Factor[((y*(l + a)~2 - 1 - a + a*b)*(x*(l + b)~2 + a*b + a) ~d - 
(x*(l + b)-2 - 1 - b + a*b)*(y*(l + a)-2 + a*b + b)-d)][[3]]; 

H = Factor [Expand [F /. 

{a -> (-1 + y + Qa)/(x-d - y) , b -> (-1 + x + Qb)/(y-d - x)}]][[3]]; 

da = Exponent [H , Qa] ; 
For [i = da, i >= 2, i — , 

H = Expand[H /. {Qa~i -> Qa~(i - 2)*(1 - x'd - y + x"d*y)}]]; 
db = Exponent [H , Qb] ; 
For [i = db, i >= 2, i — , 

H = Expand[H /. {Qb"! -> Qb-(i - 2)*(1 - y'd - x + y"d*x)}]]; 



cOO = Factor [Coefficient [Coefficient [H, Qa, 0], Qb, 0]] 

cOl = Factor [Coefficient [Coefficient [H, Qa, 0], Qb, 1]] 

clO = Factor [Coefficient [Coefficient [H, Qa, 1], Qb, 0]] 

ell = Factor [Coefficient [Coefficient [H, Qa, 1], Qb, 1]] 



uOO = Factor [Expand [cOO /. {x 

uOl = Factor [Expand [cOl /. {x 

ulO = Factor [Expand [c 10 /. {x 

ull = Factor [Expand [ell /. {x 



-> (t*y + y~d )/(t + 1)}]]; Length [uOO] 

-> (t*y + y"d )/(t + 1)}]]; Length [uOl] 

-> (t*y + y"d )/(t + 1)}]]; Length [ulO] 

-> (t*y + y~d )/(t + 1)}]]; Length[ull] 



Comment: the following code checks the positivity of the coefficients of the last factor of cqo, cqi, cio, cn 



BAD = False; 
u00[[l]]*u00[[2]]*u00[ 
For [i = 1 , i <= Length 

If[ (u00[[7]][[i]] /. 
u01[[l]]*u01[[2]]*u01[ 
For [i = 1 , i <= Length 

If[ (u01[[7]][[i]] /. 
Ul0[[l]]*ul0[[2]]*ul0[ 
For [i = 1 , i <= Length 

If[ (ul0[[7]][[i]] /. 
ull[[l]]*ull[[2]]*ull[ 
For [i = 1 , i <= Length 

If[ (ull[[7]][[i]] /. 
Print [BAD] ; 

OUTPUT: 



[3] ] *uOO [ [4] ] *uOO [ [5] ] *uOO [ [6] ] 
[u00[[7]]], i++, 

{T -> 1, y -> 1}) < 0, BAD = True]]; 
[3] ] *u01 [ [4] ] *u01 [ [5] ] *u01 [ [6] ] 
[u01[[7]]], i++, 

{T -> 1, y -> 1}) < 0, BAD = True]]; 
[3] ] *ulO [ [4] ] *ulO [ [5] ] *ulO [ [6] ] 
[ul0[[7]]], i++, 

{T -> 1, y -> 1}) < 0, BAD = True]]; 
[3] ] *ull [ [4] ] *ull [ [5] ] *ull [ [6] ] 
[ull [[7]]], i++, 

{T -> 1, y -> 1}) < 0, BAD = True]]; 
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False 



{-l + yfy\l + yf{l + t + y + y^) 



(1 + i) 



16 



A. 3 Case A = 5 

d = 4; 

F = Factor[((y*(l + a)"2 - 1 - a + a*b)*(x*(l + b)"2 + a*b + a) "d 
(x*(l + b)"2 - 1 - b + a*b)*(y*(l + a)-2 + a*b + b)-d)] [[3]] ; 



H = Factor [Expand [F /. 
{a -> (-1 + y + Qa)/(x~d 



y),b -> (-1 + X + Qb)/(y-d - x)}]][[3]]; 



da = Exponent [H , Qa] ; 
For [i = da, i >= 2, i — , 

H = Expand [H /. {Qa~i -> Qa~(i 
db = Exponent [H , Qb] ; 
For [i = db, i >= 2, i — , 

H = Expand [H /. {Qb"! -> Qb-(i 



- 2)*(1 - x~d - y + x"d*y)}]]; 



2)*(1 - y"d - X + y"d*x)}]]; 



cOO = Factor [Coefficient [Coefficient [H, Qa, 0], Qb, 0]] 

cOl = Factor [Coefficient [Coefficient [H, Qa, 0], Qb, 1]] 

clO = Factor [Coefficient [Coefficient [H, Qa, 1], Qb, 0]] 

ell = Factor [Coefficient [Coefficient [H, Qa, 1], Qb, 1]] 



uOO = Factor [Expand [cOO /. {x 

uOl = Factor [Expand [cOl /. {x 

ulO = Factor [Expand [c 10 /. {x 

ull = Factor [Expand [ell /. {x 



-> (t*y + y~d )/(t + 1)}]]; Length [uOO] 

-> (t*y + y"d )/(t + 1)}]]; Length [uOl] 

-> (t*y + y"d )/(t + 1)}]]; Length [ulO] 

-> (t*y + y"d )/(t + 1)}]]; Length[ull] 



(* proof for Delta=5, d=4 *) 

BAD = False; 

uOO [ [1] ] *uOO [ [2] ] *uOO [ [3] ] *uOO [ [4] ] *uOO [ [5] ] *uOO [ [6] ] 

For[i =1, i <= Length [uOO [ [7] ]] , i++, 

If[ (u00[[7]] [[i]] /. {T -> 1, y -> 1}) < 0, BAD = True]]; 
uOl [ [1] ] *u01 [ [2] ] *u01 [ [3] ] *u01 [ [4] ] *u01 [ [5] ] *u01 [ [6] ] 
For[i =1, i <= Length[u01[[7]]] , i++, 

If[ (u01[[7]] [[i]] /. {T -> 1, y -> 1}) < 0, BAD = True]]; 
ulO [ [1] ] *ulO [ [2] ] *ulO [ [3] ] *ulO [ [4] ] *ulO [ [5] ] *ulO [ [6] ] 
For[i = 1, i <= Length [ul0[[7]]] , i++, 

If[ (ul0[[7]] [[i]] /. {T -> 1, y -> 1}) < 0, BAD = True]]; 
ull [ [1] ] *ull [ [2] ] *ull [ [3] ] *ull [ [4] ] *ull [ [5] ] *ull [ [6] ] 
For[i =1, i <= Length[ull[[7]]] , i++, 

If[ (ull[[7]] [[i]] /. {T -> 1, y -> 1}) < 0, BAD = True]]; 
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Print [BAD] ; 
OUTPUT: 



False 



2^6 {-, , + , „ , „2 , „3\2 



(-1 + y)14/ (1 + y + y2)b ^^ ^ ^ ^ ^ ^ ^2 ^ y3) 







(l+t)35 




(- 


_l+y)13y6 


(l + y + y2)'^(l + i + y + y2 


+ .^)^ 






(l+t)34 




(- 


_l+y)13y6 


(l + y + y2)6(i + i + y + y2 


+ y')' 






(l + t)31 




(■ 


-l + y)lV 


3(l+y + y2)6(i + i + y + y2 


^ + y') 



(l+t)30 
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